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Abstract 

The method of flow equations is apphed to QED on the hght front. Requiring that the 
particle number conserving terms in the Hamiltonian are considered to be diagonal and the other 
terms off-diagonal an effective Hamiltonian is obtained which reduces the positronium problem 
to a two-particle problem, since the particle number violating contributions are eliminated. 

Using an effective electron-positron Hamiltonian, obtained in the second order in coupling, 
we analyze the positronium bound state problem analytically and numerically. The results ob- 
tained for Bohr spectrum and hyperflne splitting coincide to a high accuracy with experimental 
values. The rotational invariance, that is not manifest symmetry on the light-front, is recovered 
for positronium mass spectrum. 

Except for the longitudinal infrared divergences, that are special for the light-front gauge 
calculations, no infrared divergences appear. The ultraviolet renormalization in the second 
order in coupling constant is performed simultaneously. To preserve boost invariance we take 
into account the diagrams arising from the normal ordering of instantaneous interactions. Using 
flow equations and coupling coherence we obtain the counterterms for electron and photon 
masses, which are free from longitudinal infrared divergences. 



Abstrakt 

In dieser Arbeit wird die Methode von FluBgleichungen im Lichtkegel-Formalismus auf die 
QED angewandt. Wir konstruieren einen effektiven, block-diagonalen Hamiltonian, wobei wir 
fordern, daB diejenigen Terme, welche die Teilchenzahl erhalten, diagonal sind, wahrend alle 
anderen Terme nichtdiagonal sind. Dieser effektive Hamiltonian vereinfacht das Positronium- 
Problem auf ein Zweikorper-Problem, well die Anteile, die eine Teilchenzahlanderung verur- 
sachen, eliminiert sind. 

Mittels des effektiven Elektron-Positron-Hamiltonians, der in der elektromagnetischen Kop- 
plung von zweiter Ordnung ist, werden die Bindungszustande des Positroniums analytisch und 
numerisch untersucht. Unsere Resultate fiir das Bohr-Spektrum und die Hyperfeinaufspaltung 
stimmen sehr gut mit den experimentellen Werten iiberein. Die Rotationsinvarianz, welche auf 
dem Lichtkegel nicht mehr gewahrleistet ist, kann fiir das Massenspektrum wieder hergestellt 
werden. 

AuBer longitudinalen Infrarot-Divergenzen, die fiir die Lichtkegel-Eichung speziflsch sind, 
treten bei den Rechnungen keine Infrarot-Divergenzen auf. Die Renormierung der Ultraviolet- 
Divergenzen bis zur zweiten Ordnung in der Kopplungskonstanten erhalt man mit der Kon- 
struktion des Hamiltonians automatisch. Um die Invarianz des renormierten Hamiltonian unter 
Boosts sicherzustellen, beriicksichtigen wir auch die Diagramme, welche aus der Normalordnung 
von instantanen Wechselwirkungen entstehen. Die aus den FluBgleichungen und der Kop- 
plungskoharenz resultierenden Counterterme der Photon- und Elektronmasse erhalten dann 
keine longitudinalen Infrarot-Divergenzen. 
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Chapter 1 
Introduction 



Quantum chromodynamics (QCD) is widely accepted as the microscopic theory of strong inter- 
action, where quarks and gluons are considered as the elementary degrees of freedom. But we 
are still unable to solve this theory on the macroscopic level (in the low energy domain) and to 
obtain an accurate description of the structure of hadrons, which are the strongly interacting 
particles observed in nature. 

A central goal of modern theoretical particle physics is to build a bridge between microscopic 
(high energy) and macroscopic (low energy) domains of QCD. 

At the microscopic level the theory is defined completely, and no further development seems 
to be necessary. On the other side, when it comes to the macroscopic hadron world, there is 
no rigorous way to calculate hadron properties immediately from QCD. One still has to resort 
to different phenomenological methods and QCD-inspired models of covariant field theory to 
calculate the hadron mass spectrum, form factors, wave functions, etc. Also lattice calculations 
serving as a model-independent method in covariant field theory improved our understanding 
of the hadronic structure. However, besides different problems, this approach is still strongly 
limited by the available computer power. 

In this work we use the light-front formulation of field theory, which is best suited for solving 
relativistic bound state problems in QCD |^, because of the simplified vacuum structure. 
We believe that light-front field theory together with the renormalization group approach for 
Hamiltonians |@],[0] provides a good strategy to reach the above mentioned goal. The physical 
idea behind this approach is to use the renormalization concept for Hamiltonians on the light- 
front to get an effective, low-energy Hamiltonian. 

Recently, Glazek and Wilson suggested the renormalization scheme for Hamiltonians, called 
similarity renormalization by the authors, where they developed a renormalization group and 
the basic elements of renormalization group calculations for Hamiltonians on the light-front [Q]. 
An alternative approach for Hamiltonian renormalization, the method of flow equations, was 
proposed by Wegner 0. It is common for both methods that they renormalize the canonical 
Hamiltonian of light-front field theory to a given order of perturbation theory. A basic advantage 
of the method of flow equations in comparison to the similarity renormalization scheme is that 
one obtains an effective, renormalized Hamiltonian for a limited (truncated) Fock space. 

The method of flow equations is based on the following idea: performing a set of infinitesi- 
mal unitary transformations, with the condition that the particle number conserving terms in 
the Hamiltonian are considered to be diagonal and the other terms off-diagonal, it is possible 
to get the block-diagonal effective Hamiltonian, where the particle number in each block is con- 
served. This reduces the bound state problem to a few-body problem, since the particle number 
violating contributions are eliminated. This procedure is similar to the Tamm-Dancoff space 



truncation [33] in the sense that also in this truncation particle number changing interactions 
are eliminated. The effective Hamiltonian constructed by ffow equations is automatically renor- 
malized to the given order in the coupling constant, since an elimination of particle number 
changing sectors can not be achieved in one step but rather sequently for transition amplitudes 
from large to small energy differences. 

Low-energy QCD is challenging and explicit calculations become complicated due to the 
nonperturbative nature of this theory. To test and to illustrate our approach we consider 
QED^^i in the light-front dynamics and investigate the corresponding positronium bound 
state problem. 

This work is organized as follows: In chapter ^ we outline the theoretical framework by 
giving the key ingredients of the flow equation method and by enumerating their applications 
to problems of solid state physics and statistical mechanics done so far. The approach of 
similarity renormalization is also considered. In chapter |] we review the methods and results 
known for light-front QED^j^i. To second order in the coupling constant we obtain the effective, 
renormalized QED Hamiltonian (chapter ^, which reduces the positronium problem to a two- 
particle problem to be analyzed further analytically (chapter ^ and numerically (chapter |^) 
for positronium bound states. The renormalization issues of light-front QED are considered in 
chapter 0. 



Chapter 2 
Flow equations 



In this section we give the key ingredients of flow equations method and set up the framework 
to use flow equations for the problems of high-energy physics. 

Flow equations are introduced in order to bring Hamiltonians closer to diagonalization. 
The method is based on the numerical recipe by Jacobi, consisting of unitary transformations 
between two states which makes the connecting off-diagonal matrix elements vanishing. If this 
is repeated for all off-diagonal matrix elements again and again then the off-diagonal matrix 
elements will become arbitrarily small. It is characteristic for these equations that matrix 
elements between degenerate or almost degenerate states do not decay or decay very slowly. 
In the next section we follow the original work of Wegner to introduce flow equations for 
Hamiltonian matrices. 

2.1 Flow equations for Hamiltonian matrices 

The aim is to find continuous unitary transformation, that brings Hamiltonian to diagonal- 
ization. Continuous unitary transformation depend on the flow parameter /. We call it U{1). 
In what follows we assume that U{0) = 1 is valid and that U{oo) brings the given Hamilto- 
nian operator H or the given matrix H to the diagonal form. The operator H acquires the 
/-dependence through the unitary transformation U{1) 

H{1) = U+{l)HU{l) (2.1) 

If the transformation U{1) were known, one would find H{oo) and the problem would be 
solved. Generally, one does not know the transformation, which diagonalizes the given matrix. 
Therefore the transformation is formulated in infinitesimal form 

^ = m),m] (2.2) 

Here rj is the generator of transformation; rj is antihermitian 77+ = —77. The connection between 
the unitary transformation and its generator is given 



U{1) = exp[ / r]{l')dl' (2.3) 




where the ordering along '1-axis' is imposed. One has for the matrix elements 

dhk,q{l) 



dl p 



E(^^-,p(0 V(0 - hkAl)Vp,S)) (2-4) 



The generator t] must be chosen in a way, that the matrix H is more and more diagonal as / 
increases. We demand, that J^k^^q h'iq falls monotonously. We separate H = Hd + Hr and use 
that use that TrH"^ is invariant under the unitary transformation 

TrHl + TrHl = TrH^ = const (2.5) 

J2ky^q h\q = TrH^ monoton falls when TrH^ increases. One has 

dTrH"^ d 

~ri ~ "j7 Z^ "'q,q ~ ^ z^ '^'iA z^vlq,p'''p,q ~ '''q,pVp,q) — ^ Z^ Vp,q'''p,q['''p,p ~ '^q,q) \'^-^) 

ai, ai q q p p,q 



The right hand side must be negative. The possible choice for the generator is r]p^q = hp^q{hp^p — 
V=[Hd,Hr] (2.7) 



hg,g) or 



'^{hk,k{l) + hq^q{l) - '2,hp^p{l))hk,p{l)hp^q{l) 



We have 

dhk,q{l) 

dl p 

Jl E hlq = -JlT. hlk = -2 UKk - hq,qfhlq = -2 J] ^^ (2.8) 

k^q k k,q k,q 

Since J2kj^q h\ „ falls monotonously and is restricted from below, the derivative must vanish in 
the limit / ^ oo. Therefore one has 

r,{l) = \Hd,H]^i^^Q (2.9) 

Practically we have reached the aim, the matrix commute with its diagonal part as / ^ oo. 
The eqs. (|2.2|) and (|2.9|) are called flow equations for Hamiltonians and are the basis for the 



work presented further. 

2.2 Similarity transformation 

Another method to diagonalize Hamiltonians continuously was suggested independently by 
Glazek and Wilson [Q, which has been called similarity transformation by the authors. In this 
subsection we follow the original work of Glazek and Wilson [Q to give the key ingredients of 
this scheme. 

We define the unitary transformation, that brings the Hamiltonian operator in a form, where 
no transitions between the states with energy difference larger than A are present. The role of 
'flow parameter' plays A, which corresponds to the ultraviolet (UV)-cutoff (see chapter ^ and 
is changed continuously. 

The Hamiltonian operator H\ (and other quantities) depend then on the continuous param- 
eter A. The Hamiltonian is given as a sum H\ = Hq\ + Hix, where Hq\ is the free Hamiltonian 
and Hix contains (not renormalized) interactions and counterterms. 

We separate each matrix M into two parts, M = D{M) + R{M). Let Eix are the eigenvalues 
of Hqx, and 

_ \Ei\ — Ejx\ 
Ei\ + Ejx + A 
Uijx = u{xijx), Tijx = 1 - Uijx = r{xijx) (2.10) 



and u{x) is the function, which is 1 for small arguments and vanishes for large arguments. 
Explicit one can choose for u{x) 

u{x) = 1 for < X < - 

o 

1 2 

u{x) falls monotonously from 1 to for - < x < - 

2 
u{x) = for - <x<l (2.11) 

o 

We define D{M)ij = UijxMij and corresponding R{M)ij = rijxMij. As we introduce Uijx and 
rij\, we have the continuous transition between D{M) and R{M). It is important in order 
to avoid the divergences in renormalization equations. When one chooses u{x) = 0{xq — x), 
then one has the original definition for D{M) and R{M). The continuous transformation for 
Hamiltonian operator can be written as before in the form 

^ = te.ff.l (2^12) 

The generator r]x is chosen in a way that D{Hx) = Hx- For practical purpose to distinguish in 
the calculations between the Hamiltonian operator Hx and its 'D'-part, it is useful to introduce 
the operator Q, such that 

Hx = D{Qx) (2.13) 

Then one has for the matrix elements 

^g.,A + «^.A^ = mjx{E,x - E,x) + [vx, Hixh (2.14) 

One is not able to find Qx and rjx from this equation. The reason is, that for the given A 
the Hamiltonian operator Hx can be additionally transformed unitary without breaking the 
condition D{Hx) = Hx- Thus one has one more condition. We rewrite the above equation in 
the form 

MijA— 7^ Vij\{Ejx - Eix) = Gijx (2-15) 



where we define 



G.,x=[vx,HIxh-^Q^JX (2.16) 



Now let 



and 



«.,A^ = D{Gxh (2.17) 



and 



V^Jx{EJX-Ea) = -R{Gx)^, (2.18) 

Now all functions are defined. One has 

There are no small energy denominators in these equations. These equations should be solved 
iteratively. 

in 



2.3 Flow equations for solid state physics 

Flow equations were successfully applied to different systems in solid state physics. Unitary 
transformation in the form of exact diaginalization of Hamiltonian operator has been tested in 

(1) model of impurity in the electron band [^; 

(2) dissipative quantum systems, in particular for spin-boson problem 0; 

(3) Lipkin model 0; 

(4) problem of interacting electrons and phonons in a solid (referred to as BCS-theory) ||1CI|| . 
In all these models a system of interest couples to its environment Hamiltonian (for example, in 
(1) single impurity couples to the band of electrons; in (2) a small quantum system couples to the 
thermodynamical bath). The aim is then to decouple 'small' system from its 'large' environment 
to find a behavior of the system. It is different from what is usually done: most theoretical work 
starts off by tracing out the bath degrees of freedom and then using suitable approximation 
schemes for the time evolution of the reduced density matrix of the small quantum system. 
In general the approach of Hamiltonian diagonalization is particularly suited for studying low- 
energy properties of the system, thereby being complementary to most other approximation 
schemes. 

In some cases the aim to get the diagonal Hamiltonian operator for / -^ cxo can not be 
reached. It was shown in the original work of Wegner pf for the model of interacting fermions 
in one dimension, that a literal use of the concept of Hamiltonian diagonalization can lead to 
convergency problems. In the case discussed there the divergences appeared as / -^ oo. The 
way out, as proposed by Wegner [Q, is to bring Hamiltonian operator instead of diagonal to 
the block-diagonal form, where the number of quasiparticles is conserved in each block. 

In many cases it is enough to transform the given Hamiltonian to the block-diagonal form. 
In particular it is so, when the block, which describes the states of interest, can be treated fur- 
ther with other methods. There are many known transformations, that construct in this way 
from the initial Hamiltonian operator an effective Hamiltonian operator, acting in a smaller 
Hilbert space and which is simpler to consider. Flow equations, where block- diagonalization is 
performed, have been compared with the following transformations 

(1) the Schrieffer- Wolff transformation, which reduces the Anderson model with single magnetic 
impurity to the Kondo problem 0; 

(2) the Foldy-Wouthuysen transformation, which decouples the Dirac equation into two two- 
component equations, one of which gives in the nonrelativistic limit the known Pauli equation 

i; 

(3) the Frohlich transformation, which constructs from the electron-phonon interaction an ef- 



fective electron-electron interaction [|T0|. In all these cases flow equations re-examine the trans- 
formations used before. 

2.4 Flow equations in field theory 

In this section we set up the framework to use flow equations for the problems of high-energy 
physics. We remind, that flow equations perform the unitary transformation, which brings the 
Hamiltonian to a block-diagonal form with the number of particles (or Fock state) conserving in 
each block. In what follows we distinguish between the 'diagonal' (here Fock state conserving) 
and 'rest' (Fock state changing) sectors of the Hamiltonian. We break the Hamiltonian as 

H = Hod + Ha + Hr (2.21) 



11 



where Hq^ is the free Hamiltonian; and the indices 'd'/r' correspond to 'diagonal', 'rest' parts 
of the Haniihonian, respectively. The flow equation for the Hamiltonian eq. (|2.2| ) and the 
generator of unitary transformation eq. ( p. 7] ) are written [Q 

^ = [r], H, + Hr] + [[H,, Hr], Hod] + [[i^od, H,lHod] 

r]=[Hod,Hr] + [Hd,Hr] (2.22) 

In the basis of the eigenfunctions of the free Hamiltonian 

Hod\i>=Ei\t> (2.23) 

one obtains for the matrix-elements between the many-particle states 



dl 



[v^ Hd + H,]ij — {Ei — Ej)[Hd, Hr]ij — {Ei — Ej) H, 



rij 



Vij = {Ei - E,)H„j + [Hd, Hr]i, (2.24) 

The energy differences are given by 



n2 ni 



Ei — Ej = ^ Ei^k — ^ Ej^k (2.25) 

fc=i fc=i 

where Ei^k and Ej^k are the energies of the created and annihilated particles, respectively. The 
generator belongs to the 'rest' sector, i.e. rjij = rjrij,rjdij = 0. In what follows we use 

[dr,Hd]d = 

[dr,Hd]r^O (2.26) 

where Or is the operator from the 'rest' sector (for example Hr or fir) and Hd is the diagonal 
part of Hamiltonian. 

For the 'diagonal' (ni = 712) and 'rest' {ni 7^ 712) sectors of the Hamiltonian one has 
corresp ondingly 

—jp- = [V, Hr]dij 

„ = [V^ Hd + Hr]rij — (Ei — Ej) [Hd, Hr]rij + — (2.27) 



where we have introduced the cutoff function Uij{l) 

Uij{l) = e-(^«-^^)'' (2.28) 

The energies Ei{l) start to depend on the flow parameter / in the second order of perturbation 
theory, that is taken into account by the renormalization of single-particle electron and photon 
energies (see chapter |^. 

The main difference between these two sectors is the presence of the third term in the 'rest' 
sector J^—i22.^ which insures the band-diagonal structure for the 'rest' part 

i.e. in the 'rest' sector the matrix elements with the energy differences larger than 1/yl are 
suppressed. In the similarity renormalization scheme [^ the width of the band corresponds to 
the UV cutoff A. The connection between the two quantities is given 

I = ^ (2.30) 

19 



The matrix elements of the interactions, which change the Fock state, are strongly suppressed, if 
the energy difference exceeds A, while for the Fock state conserving part of the Hamiltonian 
the matrix elements with all energy differences are present. As the flow parameter / — >■ oo (or 
A — > 0) the 'rest' part is completely eliminated, except maybe for the matrix elements with 
i = j. One is left with the block-diagonal effective Hamiltonian. 
Generally, the flow equations are written 

— -— = [r], Hd + Hr]ij — {Ei — Ej)[Hd, Hr]ij -\ — — 

dl ' dl Uij 

,,^|J/..//J, + ^(-^^) (2.31) 

where the following condition on the cutoff function in 'diagonal' and 'rest' sectors, respectively, 
is imposed 

Udij = 1 



One recovers with this condition the flow equations eq. ( p.27|) for both sectors. Other unitary 



d\ 


— "'lA'h 




Tij 


ritj - 


E.-E, 



transformations, which bring the Hamiltonian to the block-diagonal form, with the Fock state 
conserving in each block are used 

'^ - Uij[r], Hd + Hr]ij + ^ij-TT-—^ 
and § 

where Uij + Vij = 1; and the constrain eq. ( |2.32| ) on the cutoff function in both sectors is implied. 
One can choose the sharp cutoff function Uij = 6{X — |Ajj|). 



1.-^ 



Chapter 3 

Light-front field theory 



3.1 Introduction 

The development of light-front field theory dates back to the work of Dirac |]TT|, where he 
introduced the light-front coordinates (the coordinate vector is x = {x^,x~,x±) with x^ = 
x^ ± x^ and x± = (a;i,X2)) and the concept of front form dynamics for Hamiltonians. Dirac 
suggested, that a Hamiltonian operator can 'propagate' a physical system either in the usual 
time x° (instant form dynamics) or in the light-front time x~^ (front form dynamics). The 
latter form of relativistic dynamics combines "the restricted principle of relativity with the 
Hamiltonian formulation of dynamics" [0 . 

Later the rules for front form perturbation theory were formulated |T^, and the equivalence 



of this theory with the Feynman rules of covariant perturbation theory was established [|T3[,[|l4 

Recent interest in light-front coordinates is driven mainly by two topics: low-energy bound 
state problem in QCD, where light-front coordinates offer a scenario in which a constituent 
picture of hadron structure can emerge from QCD, because of the simplified vacuum on the 



light-front [EO], pTI, |22|; and high-energy scattering processes, where light-front coordinates 



are the natural coordinates of the system |^ , |2^ . For an extensive list of light-front references 



through the early 1990's see [0, for the list of recent reviews see [^|, [^ and references within. 
Below we give briefly an introduction to light-front field theory (for introduction see also 



Light-Front (LF) quantization is very similar to canonical equal time (ET) quantization 
(here we closely follow Ref. [0). Both are Hamiltonian formulations of field theory, where 
one specifies the fields on a particular initial surface. The evolution of the fields off the initial 
surface is determined by the Lagrangian equations of motion. The main difference is the choice 
of the initial surface, x^ = for ET and x"*" = for the LF respectively. In both frameworks 
states are expanded in terms of fields (and their derivatives) on this surface. Therefore, the 
same physical state may have very different wave functions^ in the ET and LF approaches 
because fields at x° = provide a different basis for expanding a state than fields at x^ = 0. 
The reason is that the microscopic degrees of freedom — field amplitudes at a;° = versus field 
amplitudes at x"*" = — are in general quite different from each other in the two formalisms. 

From the purely theoretical point of view, various advantages of LF quantization derive from 
properties of the ten generators of the Poincare group (translations P^, rotations L and boosts 
K) W^. Those generators which leave the initial surface invariant (P and L for ET and P_, 



P±, L3 and K for LF) are "simple" in the sense that they have very simple representations in 
^By "wave function" we mean here the collection of all Fock space amplitudes. 



U 



terms of the fields (typically just sums of single particle operators). The other generators, which 
include the "Hamiltonians" {Pq, which is conjugate to x° in ET and P+, which is conjugate 
to the LF-time x"*" in LF quantization) contain interactions among the fields and are typically 
very complicated. Generators which leave the initial surface invariant are also called kinematic 
generators, while the others are called dynamic generators. Obviously it is advantageous to have 
as many of the ten generators kinematic as possible. There are seven kinematic generators on 
the LF but only six in ET quantization. 

The fact that P_, the generator of x~ translations, is kinematic (obviously it leaves x+ = 
invariant!) and positive has striking consequences for the LF vacuum[|I3|. For free fields p"^ = m? 
implies for the LF energy p+ = {m? + px) /p_. Hence positive energy excitations have positive 
p_. After the usual re- interpretation of the negative energy states this implies that p_ for a 
single particle is non-negative [which makes sense, considering that p_ = po — Ps]. P- being 
kinematic means that it is given by the sum of single particle momenta p_. Combined with 
the non-negativity of p- this implies that, even in the presence of interactions, the physical 
vacuum (ground state of the theory) differs from the Fock vacuum (no particle excitations) 
only by so-called zero-mode excitations, i.e. by excitations of modes which are independent 
of the longitudinal LF-space coordinate x~ . Due to this simplified vacuum structure, the LF- 
framework seems to be the only framework, where a constituent quark picture in a strongly 
interacting relativistic field theory has a chance to make sense ||16| , |1^, |18|, |T^ . This is the most 
attractive feature of LF-frame to approve constituent quark model desription for QCD. 

3.2 Preceding work 

As far as progress is concerned, the light-front approach is not so far along; most research effort 
has occurred since late 1980's [jl[. Progress is currently limited by conceptual issues, mainly by 
problems in renormalization program. 

Further the methods available in the light-front field theory to solve the bound state problem 
are discussed. Then the results for positronium problem on the light-front follow. 

3.2.1 Models and methods in the hght-front field theory 

The ultimate goal of the light-front field theory is to start with QCD Lagrangian and, with 
a minimum of approximation, calculate the hadron spectrum. The basic idea behind this ap- 
proach is to use Hamiltonian techniques in the coordinate system best suited for relativistic 
dynamics. For light-front field theory, physically interesting observables are quite easily calcu- 
lated then from the bound state wave-functions. 

Further we review several methods to solve the bound state problem in the light-front field 
theory. 

There are several attempts to solve for QCD bound states in the light-front approach, 
based on simplification of initial QCD^+i Lagrangian. Instead of real QCDs+i different mod- 
els, resembling QCD and having its main properties (confinement, asymptotic freedom as in 
transverse lattice model, and also chiral symmetry breaking as in coUinear model) are solved. 

CoUinear (tube) QCD model is a phenomenological model for light-front QCD, where the 
transverse momenta of all constituents are neglected. This reduces the QCD-^^i Lagrangian 
to an effective 1 -(- 1 dimensional theory, which is then solved for the spectrum, distribution 
amplitudes and form factors of mesons |2^ . 



In the transverse lattice QCD approach one formulates an effective light-front Hamiltonian 



^f, 



for SU{N) Yang Mills theory in 2 + 1 [Q ((3 + 1) |2^) dimensions using two continuous space- 
time dimensions with the remaining transverse space dimension (two transverse dimensions) 
discretized on a lattice. In the case of (2 + 1) - dimensional theory the reduction to a (1 + 1) - 
dimensional theory takes place, that enables to investigate the string tension and the glueball 
spectrum. 

The goal of such toy model studies is to build intuition which one can hopefully apply 
to QCD^^i. However, while these models have been very useful for studing nonperturbative 
renormalization in 1 + 1 dimesional LF field theories, it is not clear to what extend these results 
can be generalized to sufficiently nontrivial theories in 3 + 1 dimensions. 

In other approaches one reduces the initial Hamiltonian of light-front field theory to an 
effective, low-energy Hamiltonian, which can be solved for bound states. Different methods are 
used to construct such an effective Hamiltonian on the light-front. 

First group of methods, known as Tamm-Dancoff approach, is based on Tamm-Dancoff 



truncation [Q and projection technique in Fock space (Bloch-Fleshbach technique). By Tamm- 
Dancoff truncation one simply restricts the full Fock space to several lowest Fock sectors of 
interest. The problem of renormalization arises in this approach, since in a limited Fock space 
the diagrams, needed for renormalization but containing more Fock components in intermediate 
state as allowed, inspite of it must be thrown away. The method of iterative resolvents developed 
by Pauli, together with discretized light cone quantization (DLCQ) approach is the synthesis of 



the methods from the first group ||2^, [^. In this approach one repeats projecting high Fock 
components upon the lower in sequence several times. Thus constructed few-body effective 
Hamiltonian describes most adequately low-lying bound states. Application of this method to 
QED is discussed below. 

Second group of methods uses renormalization group concept to construct low-energy Hamil- 
tonian. Similarity renormalization scheme 0, formulated by Glazek and Wilson, develops 
a renormalization group and the basic elements of renormalization group calculations for 
Hamiltonians on the light-front. In this method continuous unitary transformation, similar- 
ity transformation, is performed to bring the Hamiltonian operator Hij to a band-diagonal 
form \Ei — Ej\ < A. It turns out, that the width of the band A corresponds to the energy scale 
and plays the role of UV cutoff, changing of which gives rise to renormalization group running 
of Hamiltonian. Also new interactions, not present in canonical theory Hamiltonian, appear 
through this scaling. 

The alternative approach for Hamiltonian renormalization, the method of flow equations [0], 
was proposed independently by Wegner. One performs a set of inflnitesimal unitary transfor- 
mations to get the block-diagonal effective Hamiltonian, with the particle number conserving 
in each block. This reduces the bound state problem to a few-body problem, since the particle 
number violating contributions are eliminated. 

There is a need for an effective formalism for gauge theories: nobody ever has solved rigor- 
ously a relativistic many-body theory in 3 + 1 dimensions. 

3.2.2 Results for light front QED^+i 

In this section we review the results obtained by others on positronium problem in the light 
front dynamics. First we discuss the flrst group of methods based on Tamm-Dancoff approach 
(Tamm-Dancoff truncation and projection technique in Fock space). The numerical recipes are 
well elaborated there. 

Brodsky, Pauli, and Tang |^ showed how to set up the positronium bound state problem 



in DLCQ (discretized light cone quantization) with a Tamm-Dancoff truncation |33] to two 

Ifi 



sectors of Fock space: (i) the electron-positron sector, and (ii) the electron-positron-photon 
sector. In order to solve the associated eigenvalue problem, the diagonalization of the discretized 
Hamiltonian, and variational methods were used. The leading results for the triplet ground 
state of positronium were obtained not quit satisfactory. To produce significant results by 
diagonalizing the Hamiltonian matrix, one would have had to include much more Fock states 
or to improve the numerical convergence of the method applied. 

Krautgartner, Pauli, and Wolz [^ derived the continuum limit in the positronium model 
discussed above for the calculation analytically, and with a prescription for handling an infrared 
divergence they showed that to leading order the binding energy for positronium results. They 
have used an effective interaction, obtained from a projecting of the |ee7 >-sector onto the 
|ee >-sector. The corresponding effective integral equation was solved numerically using the 
method of Coulomb counterterms, introduced to improve numerical convergence. The results 
obtained for Bohr spectrum and hyperfine splitting show an excellent convergence and coincide 
with the expected values. For the first time the relativistic effects as fine structure could be 
investigated (numerically) in the light-front Hamiltonian approach. 

Kaluza and Pauli |]53[ showed that reasonably accurate numerical results can be obtained 
with proper renormalization of the Hamiltonian, although logarithmic 'divergences' remain in 
their calculations. These divergences are not a serious problem if one uses a sufficiently small 
cutoff and weak coupling. For calculation of fine structure the diagonalization technique was 
improved to enlarge the Fock space feasible with the computer. But the convergence of spectra 
was rather poor, since no proper counterterms for the Coulomb singularity of the relativistic 
problem were taken. 



Kaluza and Pirner extended the work in |34| through the fine structure |36|. 



Pauli formulated the method of iterative resolvents [^ , which enables, in general, to take 
into account infinite many high Fock states by projecting in sequence high Fock sectors unto the 
lower Fock sectors. The main advantage of the method is, that it requires the inversion of the 
effective sector Hamiltonians, corresponding to the given Fock sectors, at each step of iteration, 
instead of inversion of the full Hamiltonian matrix, standing in the bound state equation. 
To avoid in the leading order the infrared divergences (coUinear singularities) the full energy 
appeared in the interaction kernel for the electron-positron bound state was replaced with a 
free energy. In the physical amplitudes the coUinear singularities are exactly cancelled by the 
dynamical terms ^^. Pauli and Trittmann have solved numerically the corresponding light- 



front integral equation for positronium bound state in the continuum limit |3^, including the 



Coulomb counterterm technique in the computer code. The results show excellent convergence 
and coincide to a high degree of accuracy with the expected values. 

Renormalization program for positronium in the light-front QED was performed using sim- 
ilarity renormalization scheme (similarity transformation and coupling coherence) by several 
authors. Perry has shown [^ that the resultant effective Hamiltonian leads to standard lowest 
order bound state results, with the Coulomb interaction emerging naturally. 

For the first time Jones, Perry, and Glazek |^9| presented using similarity renormalization 



a systematic analytic calculation in a light-front Hamiltonian approach of the singlet-triplet 
spin splitting in the ground state of positronium through order a^. The standard singlet-triplet 
splitting of positronium was obtained with degeneracy of triplet state, that recovers rotational 
symmetry that is non-manifest symmetry in the light-front field theory. 

Jones, and Perry have calculated also the Lamb shift in the same approach |JD|. The 



dominant part of the splitting between the 2Si and 2Pi energy levels in hydrogen was obtained. 
All calculations in similarity renormalization approach were done within the perturbative 
theory and results for spectrum were obtained with several approximations analytically. 
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We do not mention here applications of both methods, Tamm-Dancoff approach |iT] and 
similarity renormalization ||20|, to bound state problem for light-front QCD. 



3.3 Canonical QED Hamiltonian 

In this section we write the canonical QED^^i Hamiltonian in the light-front gauge and also 
give briefly its derivation from Lagrangian density, following the original work of Zhang and 
Harindranath [|^ . The QED Hamiltonian in the secondary quantization is given further, which 
is used in the main part of the work (in the chapters §,^,0). 



3.3.1 Canonical light-front QED'^+i Hamiltonian 

Starting with the QED Lagrangian density 



--F^^F^^ + tPii^ + e^-m)^, 



(3.1^ 



in a fixed gauge, A~^ = A^ + A^ = 0,0 the constrained degrees of freedom, A~ and ip- (V' = 
4'+ + 4'-,4'± = ^±4'] see all definitions below), are removed explicitly, producing a canonical 
QED Hamiltonian, defined through the independent physical fields A± and ip+. Details of 
derivation follow below. The resulting canonical Hamiltonian Hcan is given as a sum of the free 
Hamiltonian and the interacting term 



P- = H,an = Jdx-(fx^(no + Hi) . 
where each term in eq. ( p.2|) is written [|24 



Ho = -id' AW' A') + e ^g^ 



and 
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where {cr*} are the standard 2x2 Pauli matrices, i = 1, 2 only, e.g., a -d^ = a'd' = o'^i—-^) + 



a 



d 



), and d^ = 2d- = 2^|^. Formally one can write the last term eq. ( p.7| ) in the form 



n 



inst 



-J' 



{d^ 



rJ' 



(3.8) 



^This derivation will not include a discussion of the gauge field zero-modes. In this work we drop zero-modes. 
For a treatment that incorporates these gauge field zero-modes from the start in QED, see Esj and references 
within; (see also footnote at the end of this section). 
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where J^ = —2e^'^^. We have used the two-component representation for fermion fields intro- 
duced by Zhang and Harindranath p4| P| 

Now the details of the derivation of the canonical Hamiltonian are presented. Given C of 
Eq. ( |3.1|) , the equations of motion are 



d^F^^ = r, (3.12) 

{i^ + e4-m)tp = 0, (3.13) 

where J'^ = —eip'j'^ilj. The physical gauge A'^ = is chosen and the projection operators A^ 
and A_ are inserted into the equations of motion. Note ip^ = A^ip and ip^^^ = A^ip. Two of the 
equations are seen to be constraint equations: 

- -d+d^A- + d'd+A' = J+ , (3.14) 

id+ip. = lia'd' + ea'A' + ^^ J ip+ (3.15) 

The fact that these are constraints can be seen from the fact that no time derivatives d~ appear. 
Note a* = 7°7*. Inverting the space derivative d~^ gives 

^- = ^^" + 2^^N (3.16) 



^_ 



id+ 



ia'd' + ea'A' + ^ ) ^h 



(3.17) 



The gauge singularities in light-front QED (QCD) that arise when one tries to eliminate the 
unphysical gauge degrees of freedom by solving the constraint equations can be seen clearly 



In the two-component representation the fermion fields are given 



V'+ 







"^"^^ ^{[a'(*9^ + eA0+*m]e} ) ' (^-^^ 



where ^ is expressed through the two-component spinors Xs as fohows 



«w = i;x./^|^^%+)(Vrf""" +<.-=""') 



„t)^-i=(;) .„d x(i)^^(;) P^io) 

and s = —s. 

For completness we give the representation for the physical gauge field (transverse component) through the 
polarization vectors £(A) 

-12 „± 



^^(-) - E/^|^^KA4e--+a+,£re-^) 



9 
e(T) = ^(l,^) and e^) = -1(1, _») , (3.II) 



IP 



in inomentuin space. In momentum space, the constraint eq. ( 3.14 ) and eq. ( p.l5 ) cannot 



determine the dependent fields in terms of physical fields for the single longitudinal momentum 
k~^ = 0. In coordinate space, this implies that the A"*" = gauge has a singularity at longitudinal 
boundary. A careful treatment of the definition of l/d~^ is necessary. A suitable definition of 
(1/9+) (and consequently {l/d~^y), which determines uniquely the initial value problem at 



x'^ = 0, is given [24| 




1 /•+00 

dy'e{x' -y')f{y')+gi 



A J- 



oo 




dx 
5_e(x" - y') = 26{x' - y') , 

e(x) = 9{x) - 9{-x) . 

where f{x) is an arbitrary field with x"*" and x^ being implicitly in the argument of it. The 
functions gi, g2 and (73 are arbitrary fields independent of x~. For a discussion on these boundary 
terms see [0]. Notice that this inverse longitudinal derivative is non-local. 



In practice, we define the inverse longitudinal derivative in momentum space. We explic- 
itly put the momentum representation of the field operators into the respective terms of the 
Hamiltonian, multiply the fields out explicitly, and then replace the inverse derivatives by ap- 
propriate factors of longitudinal momentum with the restriction \p^\/'P^ > e = 0+ \P^ is the 
total longitudinal momentum of the physical state of interest]. The absolute value sign on \p^\ 
is required for the instantaneous interactions. For example, a product of two fields gives 

1 1 

■exp[— i(p — k) ■ x] — >• — —9 [\p^ — k'^\ — eV^) exp[—i{p — k) ■ x] . (3.18) 



id~^ p+ — k^ 

The dynamical degrees of freedom are A' and iIj+ ■ The canonical Hamiltonian density is defined 
in terms of these dynamical degrees of freedom 

Taking these derivatives of the Lagrangian density and combining terms, the Hamiltonian 
density takes the following simple form 



where the constraints of eq. ( p.l4| ),eq. ( p.l5| ) are assumed to be satisfied. In our 7-matrix 



representation, only two of the components of the 4-spinor ■?/'+ are nonzero. Writing these as 
the 2-spinor ^^ and inserting the constraints of eq. ( |3.16| ),eq. ( |3.17| ), Hcan takes on the form 



written earlier in Eqs. ( |3.2|) -( p77D , where surface terms such as in 

Sx^dx- ij^J^^ (tF^^) " ~ f d^x^dx'J^ (—] J+ + surface terms , (3.21) 

are dropped 0. 



■* In other words, we are defining ip^ — A+i/; = {tpi,ip2,0, 0) = {£,, 0) 

^ A set of boundary integrals (as tfie first term in eq. ( p.21| )), arising from elimination of unphysical gauge 
degrees of freedom, is associated with the light-front infrared singularity, k^ = 0. Using the above definition 
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3.3.2 QED Hamiltonian in second quantization 

In this work we use the matrix elements of canonical Hamiltonian Hcan-, given in Eqs. (|3.2| )- 
( ^■7[ ), calculated in the free basis of Hq. Below we write the canonical Hamiltonian Hcan in the 
form of second quantization. 



We use the following momentum-space representation for the field operators, [Q and [^ 
(see footnote in the previous section) 



e(x) = T.XS j '^-^^e{p^)(h,,e-^^- + dl,e^n 



^\^) = E/^|^^(4s,Ae-- + ^-c.), (3.22) 



where spinors are xi/2 = (1,0), X-1/2 = (0,1), with s = —s and polarization vectors e\ = 
^(l,z), e'^_i = 45(1, —i)', the integration running over the p+ > only these states, that are 
allowed the light-front theory. 

The corresponding (anti) commutation relations are 

[ag,A, aJ,A'] = ^g,g"^A,A' , (3.23) 

where 

<5py = 2{2nf S{p^ - p'+)8^'\p^ - p'^) . (3.24) 

The light-front vacuum has trivial structure for both boson and fermion sectors, namely a^lO >= 
0; 6p|0 >= 0, simplifying the analytical calculations. The normalization of states is according 
to 

< Pi, Si\p2, S2 >= Sp^,p2Ssi,S2 , (3.25) 

where bp,.\0 >= \p, s >. 

Making use of the field representation eq. ( |3.22[ ), we have the following Fourier transformed 
for 
the free Hamiltonian 

ff. = E / ^|^''(P")^^^^('':.^..+<.<'.J+ E/ f^''(9*)^<A<'„A . (3-26) 
the leading order 0(e)-the electron-photon coupling 

Heej = J2 d^Pid^P2d^q e [e\dq + e"dt^(hpj)p^ + h^^dtp^ + d-pj^pr + d-p^d. 



P2"'-piy 



\S1S2 ' 

xxt2^'{PuP2,-q)XsiSq^P2-p, , (3.27) 



where 



pi/ N o^ (^■P2-^^ i i(T-pi + im , . 

r (Pi,P2,g) =2— a -a ^ , (3.28) 

q+ p^ pj 



eq. (3.18) in the boundary integrals, the singularity at fc+ = is removed; the surface terms (as the second term 
in eq. ( [3.21 )) vanish from the LFQED (LFQCD) Hamiltonian and the light-front linear infrared divergences 
cancel in all physical amplitudes [p3. The problem of 'zero modes' (the singularity of A^ — gauge at /c"*" — 0) 
is hidden now in the nontrivial asymptotic behavior of the transverse (physical) gauge degrees of freedom at 
longitudinal infinity, see the first reference of |24[| . In this work we do not consider this problem. 
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We have for the instantaneous interactions of the order O(e^) 



S1S2S3S4 
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+ + ^ X 
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XS3XS4 -^ Xs\Xs2 ~ OSJS3OS2S4 + 'JS1S4OS2S3 



(3.29) 



and 



S1S2A1A2 
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(3.30) 



here 



S = S.A / + 






a_g = Ct-g.A 



^(-g+) 



6p = 6p,,^(p+), 4 = dp_,^(p+), (3.31) 

and the 5-symbol stands for the function defined in eq. ( |3.24|) , the short notation for the integral 



IS 



I^'^^I 



+ ^2^± 



dp^d p 
2{27iy 



(3.32) 
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Chapter 4 

Hamiltonian bound state problem on 
the light-front 



In this chapter we outhne the program we follow to solve the bound state problem on the 
light-front. This includes two issues: the derivation of effective Hamiltonian and solving the 
corresponding bound state equation. The first point is discussed here explicitly in application 
to positronium bound state problem. The methods to solve bound state equation (analytically 
and numerically) are discussed in the next chapters. 

4.1 Introduction 

We use the method of flow equations to construct an effective Hamiltonian starting from the 
light-front formulation which can be used to solve the bound state problem. The physical idea 
behind this approach is to use renormalization concept for Hamiltonians to get an effective, 
low-energy Hamiltonian for a limited Fock space. The bound state problem is reduced then to 
a few-body problem with the effective Hamiltonian acting on the energy scale of bound state 
formation. 

The key to renormalization of Hamiltonians is to diagonalize the Hamiltonian operator [Q]. 
We have discussed briefly in chapter 2 the methods to diagonalize Hamiltonians continuously, 
suggested by Wegner and Glazek and Wilson , which have been called flow equations and 
similarity renormalization by the authors, resp. 

It is common to both methods that they eliminate by means of a unitary transformation 
initially the off-diagonal matrix elements between states with large energy differences and con- 
tinue with states closer and closer in energy, so that off-diagonal matrix elements between states 
of energy difference larger than A are eliminated or strongly suppressed. The final aim is to 
eliminate them completely (A — » 0) and to obtain a diagonalized Hamiltonian. 

We have mentioned in chapter 3, that a literal use of this concept can lead to convergence 
problems . As was suggested by Wegner [^ , one may leave the idea of diagonalizing immedi- 
ately in favor of block-diagonalizing. If matrix-elements between states of equal particle num- 
ber are considered diagonal, then the procedure brings the Hamiltonian into a block-diagonal 
form. Application of flow equations to an ra-orbital model has shown, that procedure of block- 
diagonalizing works much better, where block-diagonalization with respect to the quasiparticle 
number (number of electrons above the Fermi edge plus number of holes below the Fermi edge) 
is performed p|. 

It becomes apparent from the calculations by Jones, Perry and Glazek [B^ on the basis of 
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the similarity transformation, that this scheme works well down to energy differences of the 
order of Rydberg, but if one goes below, then contributions in higher orders in the coupling 
become important. 

Indeed it seems to be rather difficult to obtain bound states from plane waves by continuous 
unitary transformations. In eliminating only the terms which do not conserve the number of 
particles one postpones the diagonalization, but reduces the problem to one in the space of 
fixed particle number |^. Thus for the positronium problem it is sufficient to determine the 



one- and two-particle contribution of the Hamiltonian for electrons and positrons. 

Basically the procedure is very similar to that of the elimination of the electron-phonon 
interaction |10|| which yields an effective attractive interaction between electrons responsible 



for superconductivity. Both the method of flow equations and the similarity remormalization 
||T0| (the ref. to A. Mielke) yield results different from Frohlich's original ones |Q but in very 



good agreement with more sophisticated methods. In QED it is the interaction of the electrons 
with the photons instead of the phonons which has to be eliminated. 

A basic advantage of the methods of similarity renormalization and flow equations in com- 
parison to conventional perturbation theory is, that one obtains normally less singular effective 
interactions. 

This procedure is similar to the Tamm-Dancoff Fock space truncation [33, ^] in the sense 



that also in this truncation particle number changing interactions are eliminated. 
We deflne an effective Hamiltonian 

H^ff ^ UHU+ (4.1) 

where H is the bare Hamiltonian and the unitary transformation U is determined by the flow 
equations below. 

4.2 Flow equations in the perturbative frame 

In chapter ^ we have set up the framework to use flow equations in the fleld theory. In this 
section we formulate the equations, obtained in chapter |^, for a canonical QED Hamiltonian in 
the perturbative frame. 
Flow equations read 

r,{l) = lHAl),HAl)] (4.2) 

where the Hamiltonian is given H = Hd + Hr-, with Hd and Hr including all the terms from 
the 'diagonal' and the 'rest' sectors, resp. 

Our goal is to transform the Hamiltonian into blocks with the same number of (quasi)- 
particles. This means, that we deflne the 'diagonal' part Hd as the part of the interaction 
which conserves the number of particles (electrons, positrons, photons), and the 'rest' H,. as 
the particle number changing part. In the case of QED(QCD), where the electron-photon 
(quark-gluon) coupling is present, the number of photons (gluons) is conserved in each block 
of the final effective Hamiltonian. 

As a result of the unitary transformation new interactions are induced (see below). They 
are absent at / = and are generated as / increases. They also give rise to new terms in the 
generator of transformation rjil). This in its turn generates new interactions again. 
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To be able to perform the calculations analytically we proceed in a perturbative frame 
and truncate the series assuming the coupling constant is small. In the case of QED on the 
light-front one has for any finite value of / 

H{1) = iff + ij(i) + Hf^ + iff + ... (4.3) 

where the superscript denotes the order in the bare coupling constant, H^'^'' ~ e^; the in- 
dices 'd' and 'r' indicate the diagonal and the rest parts correspondingly. The part H^ is 
the free Hamiltonian, corresponding to the single particle energies with the structure in sec- 
ondary quantization a^a, b'^b, d~^d, where a, b, d are the annihilation operators of the photons, 
electrons and positrons correspondingly; H^^' denotes the electron-photon coupling (of the 
type a^b'^b)] H^ is the second order diagonal part of the Hamiltonian, having the structure 
b'^d'^bd, b'^b'^bb, d'^d'^dd (in the light front they correspond to the canonical instantaneous (seag- 
ull) and to newly generated interactions in the diagonal sector in second order). Note, that the 
diagonal part in the flow equations is not only the free Hamiltonian but the full particle number 
conserving part of the effective Hamiltonian. The choice of only H^ as the diagonal part gives 
rise to the band-diagonal structure of the effective Hamiltonian in each 'particle number' sector 
in the similarity renormalization scheme |Q. However, this makes a difference for the diagonal 
part only if one goes beyond third order in e. 
The generator of the transformation is 

Vil) = [H,, Hr] = [Hf\ Hi'^] + [iff, if (2)] + ... = ^(1) + ^(2) + ... (4.4) 

Up to second order the flow equation reads 

^ = [,, H] = [[H;^\Hi% h!,'^] + [[H;^\ Hi% Hl^^] + mP, hI% hP] + ... (4.5) 
Also terms of higher orders in e are generated by the flow equations. 



In the basis of the eigenfunctions of the free Hamiltonian H^ 



(0) 
d 



HP\i>=Ei\i> (4.6) 

one obtains for the matrix-elements between the many-particle states 

V = iE,- E,)Hl]^ + (E, - E,)hS] + ... 

^ = -(E.- E,fH^ + [r^(^), h!:%, - (E, - E^fHi^ + ... (4.7) 

The energy differences are given by 

n2 ni 

Ei — Ej = 22 Ei,k — z2 ^h^ i^-^) 

k=l k=l 

where Ei^k and Ej^k are the energies of the created and annihilated particles, respectively. 

The energy E^ depends on the flow parameter / only in second order in the coupling. There- 
fore one has 

^EllL = ^(E.-E?H^'^- 

Hil^il) = hI^H = 0)e-(^-^^)'' = ii«(A = A -^ oo)e-^^^ (4.9) 



Here we have used the physical meaning of the flow parameter I; I = ^, where A is UV-cutoff 
(see chapter]^, eq. ( p.30| )). 

In the flow equations A defines the smooth UV-cutoff. This fact insures the analytical 
behavior of the effective Hamiltonian with A, that helps in numerical calculations. 

In second order one has to distinguish between the behavior of the 'diagonal' and the 'rest' 
term. For the 'rest' part one has 



dH. 



(2) 

"i _ r^(i) H-(i)i .. _ rj?. __ z?n2h-(2) 



dl 



[v^'\H(%.,-{E,-E,rHS, (4.10) 



where index 'r' by [q^^^ H^^']r defines the particle number changing part of the commutator. 
Introduce 

i^g(/)=e-(^-^)^'i^g(/). (4.11) 

Then the solution reads 

^;g(0 = HS^il = 0) + jyi'e^'''-''^^''[v^'\l'),H('\l')U. (4.12) 

For the 'diagonal' part one has 



dH. 



(2) 
dij 



[v^'\H^'^Uj (4.13) 



dl 
and the solution is 

hS]{1) = Hj^il = 0) + J^dl'[v^'\n,H^'\nU. (4.14) 

Note, that though in general the commutator [[H^ )H\ ], -ff^ ] is not zero, it is not present in 
the flow equation due to the definition of the diagonal part. The corresponding commutator 
[[Hfi \H^'^'^]H\ ] in the 'non-diagonal' sector insures the band-diagonal form for the 'rest' in- 
teraction and also gives rise to the different structure of the generated interaction (the integral 
term) in the 'rest' and 'diagonal' sectors. 

The commutator [ri^^\H'^'^'>] gives rise to new terms in second order in the bare coupling 
e. In the case of QED it induces new types of interactions and generates the renormalization 
group corrections to the electron (photon) masses. The coupling constant starts to run in third 
order in e. 

4.3 Effective low-energy Hamiltonian 
4.3.1 Effective electron-positron interaction 



In this section we follow mainly the work [Q. We calculate by means of flow equations the 
effective interaction between electron and positron, generated in the second order in coupling 
e by elimination of electron-photon vertex. We use the canonical QED Hamiltonian on the 
light-front in second quantization, given in chapter ^. 

The term for electron-photon vertex by finite / is given 



Asis 



1S2 



Hee^ = Y. j d'^Pid^P2d^Qi9*p,pf{^)£\ag + gp,pf{l)e"dtq){h-^J}p, + hpjtp^ + d-pj}p^ + d-pjtpj 

XX^2n(Pl,P2,-g)Xsi^g,P2-pi , (4.15) 

9fi 



and 

^i[PuP2,q) = 2— ^ '- — a -a ^ '- — , (4.16) 

g+ pj pj 

where /-dependence conies from the unitary transformation performed. Also we write exphcitly 
the momentum dependence of the couphng constant and the mass as long as / 7^ 0, here pi and 
Pf stand for the set of initial and final momenta,resp. The initial conditions for the coupling 
constant and the mass are defined at the value of bare cutoff A — * 00 (/a = 0) 

lim g{lA) = e 
lim in{liC) = m (4-17) 

A^co 

Following the procedure outlined in the previous section, the leading order generator of the 
unitary transformation is 



\s1s2 

^xt^i{Pl,P2, -q)XsrSg,p,-p, , (4.18) 

'-^PiPf "''' 
where we have introduced Ap.p^ = J2p7 ~J2 P] 5 and the light-front fermion energy is p^ = - — ^^^, 

_L2 

the photon one q~ = ^. Further we calculate the bound states of positronium. In what follows 

we consider in the |ee > sector 

the generated interaction to the first nonvanishing order 

H!i:e= E j d''pid'p2d'p^d^PiV^Xl)h^^d-;^dp,hp,xXxtx-s2XsA.+P2,P,+v. ^ (4-20) 

S1S2SS.S4, 

with the initial condition hmA-,00 K^d"(^a) = 0, 
and the instantaneous interaction 

^:l1= E /^Vc^'P2rf'P3rf'P4i^;r;;(o^+<^P2Vxixix.-.xs.5p.+p.,P3+P4, (4.21) 

SIS2S3S4, 

where 



^(0 = «(0(^^3^, 

iL-C;(^A) = ^'-' ' (4.22) 

where the initial (bare) value of instantaneous interaction (its matrix element in |ee > sector 
eq. ([4.22|) ) is defined by the instantaneous term T^geee ^q- ( p.7|) ,eq. ( p.8|) of the canonical light- 
front QED Hamiltonian. The order of the field operators in both interactions satisfies the 
prescription of standard Feynmann rules in the |ee > sector. 

We note, that generated interaction eq. ( [4.20| ) is a new interaction, induced by flow equations 
in the second order in coupling, and corresponds to the dynamical photon exchange; while the 
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instantaneous term eq. ([4 .211) enters the canonical light-front QED Hamiltonian and describes 
the instant photon exchange. Instantaneous term is special for the light-front calculations Q. 

To the leading (second) order we neglect the / dependence of the energies (which start to 
run in the second order in coupling constant, see chapter |^) in the interactions, that enables to 
write the flow equations for the corresponding couplings. 

The flow equations to the first (for the electron-photon coupling) and second (for the in- 
stantaneous and generated interactions) orders are 



dl P^Pf 

dl 



^PiPfdpiPfK'') 

(4.23) 



where 

\.Pf=T.Pi-J2Pj (4-24) 

The second and the third equations are written for the 'diagonal' sector. The matrix ele- 
ment < [r]^^\l), Hee-y] >\ee> IS uudcrstood as the corresponding commutator between the free 
electron-positron states, namely < P3S3,p4S4^\...\piSi,p2S2 >■ In the order O(e^) there is also 
the contribution of the commutator < [q^^^ , -^667] > to the free Hamiltonian in one-electron and 
one-photon sectors, that defines the renormalization of electron and photon masses, resp., (see 
chapter [Zp. Renormalization group running of both (instantaneous and generated) interactions 
starts in the order O(e^), and the electron-photon coupling starts to run in the order 0{e^). 
Neglecting the dependence of the light-front energies on the flow parameter / (that is the 



corrections of higher orders), the solution of eq. ( |4.24| ) reads 
^P,P/(0 = fp^Pf■e + 0{e^) 
d(0 = g'pXilA = 0) = e' + 0{e') 

KTfi^) = ///'<[^^'nn,^ee,(0]>|ee->+O(e^) 

_a2 ; PiPf 

fp,Pt = e-^-'-f' = e-^^ , (4.25) 

where the subscript |ee > means, that the commutator is considered in the electron-positron 
sector. The electron-photon interaction exists in the band of size A (|Ap^p | < A), whereas the 
matrix elements of instantaneous and generated interactions in |ee > sector are defined for all 
energy differences. 

We give below the explicit expressions for the generated interaction, and details of calcula- 
tions can be found in Appendix 0. In what follows we use the notations of this Appendix. 

The matrix elements of the commutator [r]^^\ H^e-y] in the exchange and annihilation chan- 
nels are 

2m fp+_p+\ ylpi,P'J,9p4„P2 ~r Vp4,P2 9pi,P3) 1 

<W'\H,,^]>={ ' ' (4.26) 



r(an) 1 



^•^211 (p+_|_p+) WPl.~P25'p4,-P3 ' Vp4,-P39pi, 



-P2) 



^ The analogous problem, interacting electrons and phonons in a solid, was considered by Wegner and Lenz 
pO| . They used 'equal time' canonical electron-phonon Hamiltonian known from BCS-theory. The elimination 
of electron-phonon interaction by the flow equations generates a new interaction between electrons, that defines 
an effective electron-electron interaction. This effective interaction leads to the leading Coulomb behavior. 

9,8 



where 
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dl 




9pi,P2\'') 
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JP1,P2 


(0 






^Pl,P2 ~ 


Pi 
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-(Pl 
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P2, 



(4.27) 

and the conservation of "+" and "_L" components of the total momentum is imphed, i.e. 
Pi +pt = pt ^Pt ^^d Pi +P2 = Ps" +P4 ■ The matrix elements M2M between the corresponding 
spinors in both channels are given 

^2? = [xtPiPuP^iPl - P3)XsA [Xt2^\-Pi, -P2, -{Pl -P3))Xsi] 

(4.28) 

^iT^ = [xt^\-P4,P3,-{Pl+P2))Xrsi][xt2'^\Pl,-P2,Pl+P2)Xsi] 

where 

r(PuP.,)=2^^-'-^i^a'-a"LJ^ (4.29) 

q+ pj pj 



This equation, eq. ( [4.29|) , defines the spin structure of the generated interaction. 



We combine the formulas for commutator [ri^^\ H^e-y] together with the generator t]{1) and 
coupling constant g{l), expressed through the similarity function /(/). The generated interac- 
tions eq. ( 4.25 ) are given then in both channels 

\Pl PS ) \ ^Pl,P3 ^P4,P2 



-I- / JA d\' JP4-P3,X'"'^ Jx dX' Jpi-P2,A'"A 

[Pl +^2)1 ^Pl-P2 ^P4-P3 



Yi<^'^)(X) = e^M^''"^ ^ ' ■^'^ "'^^'" •^P4-P3,A'"^ I Ja "dA 



(4.30) 



where in the integral we have neglected the dependence of light-front energies on the cutoff 
A (that is the correction of the order 0(e^)), and the connection between flow parameter and 
cutoff, / = 1/A^, is used. We use the explicit form for similarity function 



P1P2 



fp.,P2,x = e"^^ . (4.31) 

that gives for the generated interaction in both channels fig. (2), fig. (3) 

T/(e^)('\'l — ^2/i^(ea:) 1 ^Pi,P3 + ^P4,P2 /i _ / f \ 

^ qen l^J ~ ^ ^'^2n / + +\ a 9 , A9 V J pi;pz,Xj pi,p2,X) 



,2 /i^(ea:) ^ ^Pi,P3 "T ^P 4,P2 

{pt-pt) Al^^.^ + Al^^^ 



(4.32) 



\/(an)(^\ _ 2 »^(an) ^ ^Pl,-P2 + ^P4,-P3 /i _ f f \ 

^gen \^J ^ ^^■^2ii /„+ i „+\ a2 i a2 V J pi,-p2,XJ p4,-p3A) ' 

\Pl ~r P2 ) ^Pl,-P2 "T ^P4 -P3 



For finite values of / the solutions of the flow equations have no divergences in the form of 
small energy denominators, present in perturbation theory. The divergent contribution in the 
generated interaction as Ap^p^ ~ '^P4,P2 ~ is effectively cancelled by the factor in bracket 



9.P 



containing similarity functions (1 — /pi,p3,a/p4,p2,a) (the same for annihilation channel). One has 
for the generated interaction in the exchange channel 

^Pl,P3 "I" ^P4,P2 _ ^P1,P3 ' ^P4,P2 /A oo\ 

A2 I A 2 ~ A2 I 2A A 

Pl,P3 ^ ^P4,P2 '-^PiPf ^ '^^Pl,P3^P4,P2 

where 

^PiPf = Pi + P2 ~ P3 ~ Pa — ^pi,P3 ~ ^P4,P2 — ^pi,-P2 ~ ^P4,-P3 (4-<j4) 

due to the total momentum conservation in '+' and 'transversal' directions. The matrix ele- 
ments with any energy differences (i.e. VAp.p^) are present in 'diagonal' sector. 

The effective Hamiltonian is defined in the limit A — »• 0. In this limit the electron-photon 
coupling, present in generated interaction through the similarity functions fp^pf, is completely 
eliminated fp^pfiX — > 0) = for Ap.p^ 7^ 0, and generated interaction is given by the expression 
that does not depend explicitly on the cutoff A. 

The resultant generated interactions in both channels are given fig. (2), fig. (3) 



gen .,. ^2 ^ ^2 

y,en - e iV2,, ^, ^ ^,, , (4.35) 

where we have introduced 

Ai _ Ai . A _ A2 



A _ Ai_ _ Ai . A _ A2 



^PlP3 P+ {x'-x)P+ ' ^P4P2 P+ (x'-x)P+^ 

^ - Ml . ^ - Ml f4 361 

^Pl,~P2 ~ P+ 1 ^P4,-P3 ~ P+ {^.OV) 

(see Appendix ^ for the explicit definition of these quantities in the light-front frame). 

The expression eq. ( [4.35|) is written for the rescaled value of the potential, i.e. Vx = Vx/P^"^, 
and the cutoff is defined in units of the total momentum P^, i.e. A -^ p+ , with / = 1/A^. The 
spin structure of the interaction is carried by the matrix elements M2ii, defined in Appendix ^ 

We summarize the instantaneous interactions in both channels in the order O(e^) 






fig. (2), fig. (3) 

~{p\-ptY """""" 

{an) 4e2 

^inst ~ / + I +\2 "siS2"S3S4 5 V^-'-'') 

[Pl +P2) 

where we have used xt3Xt2^XsiXs4, = '^siS3^s2S4 + '^sis2^s3S4- For the rescaled potential in the 
light-front frame Appendix ^ eq. ( [A.14| ) we have 



T/.('=^) - zr A A 

^mst / _ A2 "siS3"s2S4 



V^ = Ae^ 6s,-sA.-s. ■ (4.38) 



."^n 



In the second order in coupling the effective electron-positron interaction, calculated in the light- 
front dynamics, is given as a sum of interactions, generated by the flow equations eq. (|4.35| ), and 



instantaneous interactions eq. ( [4.38| ), that are present in the light-front gauge calculations;i.e. 

V^ff = V^^- + V^-^'= Y. iy^l + v!:nst) (4.39) 

(i)=ex,ann 

where we sum over exchange and annihilation channels. 

In the chapters § and ^ we use the effective electron-positron interaction eq. (|4.39| ) to 
calculate positronium mass spectrum. 

In the next section we outline the general computational strategy to solve for bound states 
of positronium using an effective Hamiltonian. 

4.3.2 Positronium model (general computational strategy) 

In this section we give the program to solve positronium bound state problem in the light-front 
dynamics. This approach can be applied also to the other systems. 

Required that the particle number (Fock state) conserving terms in the Hamiltonian were 
considered to be diagonal and the other terms off-diagonal an effective Hamiltonian was ob- 
tained which is block-diagonal in particle number (Fock) space. This means, that the 'diagonal' 
sectors are decoupled, since the particle number (Fock state) violating contributions are elim- 
inated, and one is able then to truncate the full Fock space in the effective Hamiltonian to 
the lowest Fock sectors of interest. The bound state problem is reduced then to a few-particle 
Hamiltonian problem. Sure, this is true to the given order of perturbation theory. This is the 
idea of the approach. 

Below we illustrate schematically this procedure for the case of light-front QED, where the 
positronium bound state problem is reduced to a two-particle problem. 

We start with the light front Schrodinger equation for the positronium model 

i/iclV'n >= M„Vn > (4.40) 

where Hlc = P'^Pfi is the invariant mass (squared) operator, referred for convenience to as 
the light front Hamiltonian of positronium and \ipn > being the corresponding eigenfunction; 
n labels all the quantum numbers of the state. 

The canonical Hamiltonian of the system H^c contains infinitely many Fock sectors (i.e. 
one has for the positronium wave function \ipn >= Cee|(ee)„ > +Cee-y\{ee'y)n > +Cee77l(ee77)„ > 
-|-...) and each Fock sector contains states with arbitrarily large energies. We now 

(1) introduce the bare cutoff (regularization) with the result H^{A) - the bare Hamiltonian; 

(2) perform the unitary transformation by means of flow equations with the result if ^^^-^ - 
the effective renormalized Hamiltonian (table 1 for finite value of A); 

(3) truncate the Fock space to the lowest Fock sector (|ee >) with the result H'^-^-f^ - the 
effective, renormalized Hamiltonian acting in the electron-positron sector. 

Then the eigenvalue equation reads 

if^^^|(ee)„>=M2|(ee)„>. (4.41) 

The effective light front Hamiltonian consists of the free (noninteracting) part and the effective 
electron-positron interaction 

ffeff ^ ^(0) ^ yeff ^4^2) 



The light front equation eq. (^4.41 ) is then expressed by the integral equation (the coordinates 
are given in fig. (4) 

w? + k'l 2 1 / / "' 

^,il_^,\ - M^ I Mx , k^] S3, S4) 

r (JTrj Ic I — */ — * — * 

+ J2 0/0 A^ < x',kj_;s3,S4\V^^\x,k±;si,S2> ipn{x,k±;si,S2) = 

(4.43) 

The kernel of this equation, the effective electron-positron interaction V'^-^'^ , was obtained in 
the previous section, eq. ( [4.39| ). Note, that the effective interaction V^^^ , present in eq. ( f4.43|) 



is boost invariant (i.e. does not depend on P^). The integration domain D in eq. ( [4. 431) is 



restricted by the covariant cutoff condition |^ 

771 -\- h" 

— ^ < A^ + Arn^ (4.44) 

x[l — x) 

which allows for states which have a kinetic energy below the cutoff A. 

In the chapter || we solve eigenvalue equation eq. ([4.41|) analytically, using the bound state 
perturbation theory. In the chapter ^ we solve the corresponding integral equation eq. ( |4.43| ) 



numerically, using the numerical methods elaborated in Tamm-Dancoff approach and discussed 
in the section 13.2.2. 



We have considered here |ee >-sector. In the next section we consider extended Fock space 
to give diagrammatic representation of the effective Hamiltonian, obtained in the second order 
in e. 

4.3.3 Effective, renormalized QED Hamiltonian 

In this section we review diagrammatic rules for the effective Hamiltonian, obtained by flow 
equations in the second order in e. It is useful to represent Hamiltonian matrix elements in 
the form of table, consisting of different Fock blocks, each of them has inflnite many matrix 
elements describing the transition between corresponding Fock states and with inflnite many 
energy differences (Table 1). In this representation the canonical fleld theory Hamiltonian on 
the light-front has deflnite structure, that is simpler than that in equal-time formalism. On 
the light-front the vacuum 'does not fluctuate'. This means, that only one- and two-particle 
transitions are possible, three- and more-particle transitions are absent. In the case of QED 
on the light-front one-particle transitions are deflned by the electron-photon interaction and 
two-particle transitions are given by the instantaneous interactions. Kinetic energy term is 
diagonal in Fock space representation, also there are instantaneous terms that does not change 
particle number. The canonical QED Hamiltonian has therefore pentadiagonal structure in this 



representation [2S 



Now consider the effective QED Hamiltonian. We do not include instantaneous diagrams, 
since the flow equations do not change them to the second order in coupling. We start with the 
situation, when the 'rest' (Fock state changing) sector consists of matrix elements of electron- 
photon vertex, and the 'diagonal' (Fock state conserving) sector has matrix elements of kinetic 
terms for electron and photon. We perform the unitary transformation to eliminate the 'rest' 
sector. In the second order in e the elimination of electron-photon vertex generates one- and 
two-particle operators. The elimination of 'rest' sector in the next orders in e generates many- 
particle operators. 



The matrix elements of the effective Hamiltonian, obtained by flow equations in the second 
order in e, namely the diagrams in different Fock sectors are depicted in Table 1. Dot denote 
the zero in the second order matrix elements. Corresponding analytic expressions for the 
matrix elements in 'diagonal' and 'rest' sectors are listed in fig. (2). The diagrammatic rules 
are obtained by direct calculation of matrix elements between free particle states. We consider 
the situation with finite A. The matrix elements of the 'rest' sectors are squeezed in the energy 
band A^-p = \J2p7 ~ J2p~f\ < -^S there are matrix elements with all energy differences in 
'diagonal' sector. As A — >■ the 'rest' sector is completely eliminated to the given order of 
perturbation theory. One ends up with the block-diagonal effective Hamiltonian, where in each 
block the Fock state is conserved. 

For completenes we give below the matrix elements of effective interactions in different 
Fock sectors, generated by the flow equations in the second order in coupling constant. The 
transitions between Fock states and corresponding matrix elements are given, resp., 
in the 'diagonal' sectors 



lee >^ lee >, \eeee >- 
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\eeee >, ... 

poo d-fpip-^x' f i\i coo 4fp4P2A' f i\;' 
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pipiiX 
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P4P2A 



ee7 >— > |ee7 >, |ee77 >— > |ee77 >, ... 
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in the 'rest' sectors 

|ee >^ \eeee >, \eeee >^ 



^P2k2X 



ee >, 



(4.45) 
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, (4.46) 
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where 'dots' denote the higher Fock sectors; one obtains the next higher Fock sector, when in the 
given Fock state an additional (ee)-pair or 7-photon are created. Note, that in eq. ( [4.45|) ,eq. ([4.46|) 
the order of momenta is given for |ee > channel, fig. (2). We use the similarity function 



pjpf 



fp.PfX = e A2 



(4.47) 



in eq. ( |4.45| ),eq. ( [4.46| ) to get the explicit form of interactions in both sectors listed in fig. (2). 
'Rest' diagrams are drawn schematically to show the difference between the interactions in 
'diagonal' and 'rest' sectors. Namely, for the 'rest' sectors we imply, that the corresponding 
momentum exchange must be done in the diagrams of fig. (2) to get analytical expressions for 
diagrams depicted in Table 1. 

In the second order flow equations give rise to the mass corrections in one-particle sector. 
It turns out that the electron, photon mass corrections are equal, but have the opposite sign. 



.-^."^ 



to the standard electron, photon self energy terms, obtained in the light-front perturbation 
theory. We consider this explicitly in the chapter |^. In order to make the diagrammatic rules 
complete we include mass diagram for the photon in the table (sector I7 >— i> I7 >). After the 
unitary transformation is performed the electron (photon) mass depends on the cutoff 

ml = ml- 5Sa , (4.48) 

where 5Tj\ is the self energy term (mo = for a photon), and subscript denotes the bare mass. 
In the third order the coupling constant gets the cutoff dependence, i.e. e\ = eo(l + 0(e|)), 
that is beyond our consideration. 

The two-component LF theory, introduced by Zhang and Harindranath p4|, as compared 



to four- component formalism of Brodsky and Lepage, is formulated purely in terms of physical 
degrees of freedom; so that each term in the effective, renormalized Hamiltonian corresponds 
to a real dynamical process (or give renormalization term). Different 'diagonal' sectors of the 
effective Hamiltonian contribute to: in I7 > sector - to self energy photon operator, in |ee > 
sector - to electron-positron bound state (or corresponding scattering process), in I77 > sector 
- to light-light scattering, in |ee7 > sector - to Compton scattering, et. cetera (see Table 1). 

In the previous section we considered this formalism in application to the positronium bound 
state problem. 
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Chapter 5 

Positronium spectrum (analytically) 



In this chapter we solve an effective eigenvalue equation for positronium, obtained in the pre- 
vious chapter, analytically. We perform the bound state calculations perturbatively. The idea 
of calculations is the following. We split an effective Hamiltonian H'^-^^ into H^^\ a part which 
is solved nonperturbatively, and 5V, the difference between the original Hamiltonian and H^^\ 
The effects of 6V are to be computed using bound state perturbation theory. The criteria for 
choosing H^^^ is that it approximates the physics relevant for the given bound states (positron- 
ium in our case) as closely as possible. As a consequence, H^'^^ contributes the dominant term to 
the mass spectrum and the bound state perturbation theory converges with respect to SV/H^'^\ 
For QED, where the analytic answer is known, the leading order solution H^'^^ is simply given 
by a sum of kinetic terms and the Coulomb potential. This result arise straightforward from 
the form of the effective electron-positron Hamiltonian in the nonrelativistic limit. For more 
complicated theories as QCD the hint to choose H^^^ comes from phenomenological models. 

In the next section we define explicitly bound state perturbation theory for positronium 
system. The calculation of Bohr spectrum and the ground state spin-splittings is given further. 



In analytical calculations of singlet-triplet spin-splitting we follow the work |39 



5.1 Bound state perturbative theory (BSPT) 

First introduce instead of the light front parameterization fig. (3), used before for the single- 
particle momenta, the instant form. We express the variable {x, k±) in terms of the equal-time 
variable p= {pz,K,±) as 

f=pl + Hl (5.2) 

and similarly for x' and p ^ as function of p^. 
The Jacobian of this transformation, J{p), is 

^^^^ ^Wz^ 2(/+m2)3/2 • (5-3) 

The instant form is used for practical purposes: it is simpler to recover the rotational 
symmetry there, the symmetry that is not manifest in the light-front frame. 
We choose the leading order Hamiltonian operator for positronium 

i7(0) = h + v^^^^ ^ (5.4) 



where h is the free part (sum of corresponding kinetic terms) and Vcoui is the Coulomb interac- 
tion 

Let us solve the corresponding Schrodinger equation on the light-front 

i7(°)|V'^(P)>=^^|^jv(P)>, (5.6) 

where P is the positronium momentum. The eigenvalues and eigenf unctions for positronium 
bound state on the light-front are defined in a standard way 

^N — — p+ — 



IMP) >= Es,s, Ip,p, VP^P^2(27r)3 5(3)(P-pi -p2)$iv(x«:xSiS2)6+(piX(p2)|0 > 

^^^1^2 2{2jy" "" ^*Nix'^±siS2) ^'n{xk±siS2) = 6nN' (5.7) 

Mtv stands here for the leading order mass of positronium. Combining the definitions for the 
wave function and the energy together with the light-front Schrodinger equation, one has 
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or, after change of coordinates according to eq. 



/ o -2 , \ -. ^ I d^P\/Jip)Jip') 

[Ml - A{p' + m^)] ^n{p's,s,) = Y: 2(27r)3 ^-"'(^^'^O "^^(^^1^2) , (5.9) 

where the wave functions are redefined to have the norm 

J2 [ d^P^*N{psiS2) $'iv(psiS2) = 5nn' ■ (5.10) 

S1S2 •' 

Our aim is to obtain the nonrelativistic Schrodinger equation for positronium. Note, that in 
the nonrelativistic limit jp' /m? « 1 we have 

/WW « i (1 - i^ilf^) 

Mn = {2m + BNf ~ Am^ + AmB^^^ , (5.11) 

where the leading order binding energy B^ is introduced. Then in the leading order the bound 
state equation for positronium is 

f ^2 
pi 

m 

(5.12) 



-BAi,(iys,s,) = -y:j,^,(J--±^i-v.^;^ *„(p.,.,) 



Using the explicit form for the Coulomb potential, eq. ( ^.5|) , we get the equation that 
determines the leading order bound state wave function: 

--5J ^,{P') = A / 7^-^,^,iP) (5-13) 

m I 277^ J [-p — p'Y 



with 

^N = $M,Se,.e-(P'S3S4) = <^^,{p')Ss,s3Ss,s,- (5.14) 

This is the standard nonrelativistic Schrodinger equation for positronium. The solution is 
characterized by ^ = {n,l,m), the usual principal and angular momentum quantum numbers. 
The wave functions are given through the hyperspherical harmonics 

'^ ^ Aef ' 

Yfx = yn,l,m = fn,l{^)yi,mid,(l)) 

ma^ ma 

B^ = -^, e. = — (5.15) 

and for the binding energy one has the standard nonrelativistic expression for positronium 
bound state to O(e^). We write for completeness the coordinates used in the solution 

{el = -mBN,p) — > iuo,u) 

Mo = cosu; = 

e„ + P^ 

P ■ 2e„p 

u = -p^smt^ = , (5.16) 

\P\ ei + p^ 



but, for details, refer to 

Now we define BSPT to solve positronium bound state problem. We introduce the potential, 
arising in the nonrelativistic Schrodinger equation, eq. (|5.13| ), 

V\pfs2,Si]psiS2) = lim — - — V^" 



K«i 



2(27r)3 Am 
1 1 



2(27r)3 2m4m J? 



lim (Kxch + Knn) . (5.17) 



«i 



where V^q = Vexch + Vann is the effective electron-positron interaction, obtained in the previous 
chapter. The leading order solution is given in eq. ( [5.151 ). We perform perturbative bound state 
calculations with respect to the difference 

a 1 

SV = V'{p's3S4;psiS2) - {-TT^)— — ^77'^3i33^^2.4 , (5.18) 

27r^ [p - p')^ 

Note, that, in order to define the Coulomb potential, i.e. the ee interaction in the leading order 
of BSPT, we take only the first term of nonrelativistic expansion of the Jacobian J{p). 
Further we use the matrix elements of 6V, that are defined as 

< <^nlm\SV\^nlm >= J d'pd^p'<^li^{p)6V<!>nlm{p') , (5.19) 

where ^nim are the Coulomb wave functions. 
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5.2 Effective electron-positron interaction in light-front 
and instant frames 

We summarize together all the terms defining the electron-positron interaction V^^^ , eq. (|4.39| ) 



V'ff = Ve.ch + Vann= E {V'''' + V 



inst\ 



(5.20) 



channel 



we remind that tilde above the interaction terms denotes the rescaled potential {V = V/P~^'^), 
that does not depend on the total momentum P~^, i.e. is invariant under the light-front boosts. 
In the light-front frame the generated interaction and instantaneous term are given resp. 
in the exchange channel 



ygen 

■Trinst 



-e^Ni 



Ae' 



O.Ci Co Oq 



(a; _ a;/)2 "^1^3 -S,S, 

in the annihilation channel 



(5.21^ 



-\rinst 






M^ + Mo^ 

^6 Osis^ ^S2S4 ) 



(5.22) 



The functions Ni,N2 (current- current terms) and the energy denominators Aj,i = 1,2,3, 
Mq , Mq^, M^ are defined in the light-front dynamics ^^ as follows (see Appendix ^ fig-(3) 



-'''1 — "S1S3OS2S4-' 1 ■-'2 ^SiS2^SiS3^S2S4,'^'^ 



[X — X 



/\2 



xx'{l — x){l — x') 



+imy2{x' — 



X) 



X _ X -^1 rp± ± _|_ r r -^ rp± 

0siss0s2S4 ,-'1 ■ ^si "T "S1S3OS2S4 /^ x/^ ,x-'2 '^ 

XX 1 — X 1 — X 



^2 



iVc 



2 — '^siS2"s3S4-'3 ■-'4 ~r 0siS2^S3S4,0siS3'^^^ ,/^ n/^ 

XX 1 — X 1 



(l-x)(l-x') 



+imv2 



Sl 



S3 



OS3S4OS1S2 1^ X-'S ■ ^Si '^S3S4"SlS2 //T /\ 4 '^54 

X 1 — X X 1 — X 



and 



Ti^ 



T,* 



71- 



^ {k^-k!^)' ^ k\{s2) ^ k1(s2) 



(5.23) 



(x — x') (1 — x) (1 — x') 



^ (x — x') X x' 



^l\S?^ . '^±{S3) 



i^lisi) i^±{si) 



(1-x') ' " (1-x) X 



"' X' (1-x')' ^ 

/€_l_(^Sj = Kj_ + lS6ijKj_ ] Eij = £ij3 

with the definitions 

{xk'j_ — x'k^Y + m^(x — x'Y 



s = —s 



Ai 



XX' 



A2 = AiL_ 



>(l-x),x'^{l-x') 
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x' — X 


; A 

m2 


x(l- 


X)' 



x'{l — x') 
p-^y^^__ZK. p=(^p+^p±y^ MN = 2m + Bn. (5.24) 

Here x is the light front fraction of the electron momentum, P is the total momentum of 
positronium and B^ the binding energy of the positronium. The effective interaction eq. ( |5.2CI|) , 
generated by the flow equations, is defined in the whole parameter region, (except maybe for the 
Coulomb singularity point q = p—p' = 0, where we are not able to eliminate the electron-photon 
vertex) as follows 

2„ /A, +A.A / ie^ , , \ 






(5.25) 



To get the effective interaction in the instant frame the substitutions x{pz),x'{p^) eq. (|5T 
are to be done, also in the instant frame holds 



M^ = A{]f + m') ; M^^ = 4(p' + m^) ; (5.26) 

At the end of this section we illustrate, that the effective interaction eq. (|5.25|) gives in the 



leading order of nonrelativistic expansion \p\/m <^ 1 the Coulomb interaction. In the leading 
order one has (for the exchange channel, that gives the dominant contribution) 

Ai ~ A2 = A = {p-p'f 



V^ 



en^ . Nl 



{p-p'y 



A = ^1^^ (5.27) 

x' — X 

and the electron-positron interaction is 

yl-> ~ -e^ \ - , ^"^ ,J s,sJs,s, (5.28) 

Using the following expressions 






\2 



yx - x'y 

{p- p'f = {k^_ - k'^Y + {p, - p'^y ^ {k^_ - /?!)' + 4m^{x - x'y (5.29) 

one obtains in leading order of the nonrelativistic approximation the 3-dimensional Coulomb 
interaction (e^ = Aira) 

V\''> ^ 16m' (- ^^ ] Ss.sA.s, (5.30) 



Hence the rotational invariance is restored in this order. This result eq. ( ^.30|) does not depend 
on the details of the similarity function /a (A). 
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5.3 Positronium's ground state spin splitting 



In this section we use the effective electron-positron interaction eq. ( ^.25|) to calculate the ground 
state singlet-triplet splitting for positronium. We follow the work [^, where the similarity 
renormalization scheme was used to get an effective electron-positron interaction. 

There is an important difference between the two approaches, flow equations and similarity 
transformation. Flow equations eliminate the matrix elements between states with large energy 
differences {\Ei — Ej\ > A, where A is an effective UV-cutoff), but only for those blocks that 
change the number of quasiparticles. The value of A = A ^ cxd corresponds to the initial 
bare Hamiltonian, a flnite A determines the effective Hamiltonian at an intermediate stage, for 
A = the elimination of the particle number changing sectors is complete and the effective 
Hamiltonian is block-diagonal in particle number. One can work then in a few (or even one) 
lowest Fock sectors. 

The effective Hamiltonian, obtained by the similarity transformation, is band-diagonal in the 
'energy space'. ^ The width of the band A (namely A^/P+) introduces the artificial parameter 
in the procedure, which must be adjusted from the physical reasoning. The effective UV-cutoff 
must be low enough to neglect the contribution of high Fock states, but is restricted from below 
to stay in perturbation theory region. In the case of positronium the "window of opportunity" 

^, that is 



p+ 



and 



is simple to choose, since there are two dynamical energy scales in QED 

one of the reasons why QED calculations have been always so successful. Namely one chooses 



m^a^ 



<< A^ << m^a, and the effective electron-positron interaction does not depend on A and 



is defined on the energy scale of positronium bound state formation |]5D|. In the case of QCD, 
unfortunately, there is no such 'window', and the procedure of fitting the cutoff seems to be 
nontrivial |^D|]. In the case of flow equations elimination of particle number changing sectors 
can be performed completely, so that there is no artiflcial parameter (as effective UV-cutoff) 
left in the effective Hamiltonian. 

We use the potential eq. (|5.17|) , that appears in the nonrelativistic Schrodinger equation for 
positronium eq. ( ^.13| ), and deflnes the nonrelativistic binding energy B^ {^n ~ 4m^-|-4r?7,i?Ar). 
The effective electron-positron interaction eq. ( |5.25 ) has the following form in the nonrelativistic 
limit 



v 



1 1 






2(27r)3 4m2mV 2m2 

^ '^exch ' ^ ann 



V-ff 



{p-p'Y {x-x') 



,\2 "SiS3 0s2S4 



+ i;^ + ^'' ^'^'-^ ^' 



2^4 5 



(5.31) 



where the energy denominators were simplifled as 






2m 



1 + 



2711? 



Ai = As = {p- p'f + 0\m 



Ai = As 



2m{p' — pY 
(P'z - Vz) 




A 



2m{p' — p)"" 
(P'z - Vz) 



^Under the 'energy space' we understand the basis of the free HamUtonian Hq. 
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r2 _ !\/r'2 _ A^2 , n, [ ^2 I P 



M' = M' =4m' + 0im' 



IQ — IVlQ 



m 



The expression for Jacobian of the coordinate change is given 



J{p)J{p') 



1 
2m 



2mP' \ m? ' m? 



Combining all together, we have 

^^'^ ' 2(27r)3 4m2m V 2m^ ) 



V^ ^5v7 7 TV2 \ "' 2 J ■^I'^S "S2S4 



+ 



Am? 



+ 4e=^ (5,,52 (5, 



S3S4 



(5.32) 



(5.33) 



(5.34) 



We expand the factors A^i and N2-, appearing in the interaction, in the nonrelativistic limit. 
The term A^i contributes in V' in the order 



TtTi 



16m2% ( 1 + ^ ) + 16^ [n\p, + kIp 



Qz 



m,^ 



Qz 



z 



' \2 



-16i(si + S2)[k,^k,±] - 4(/t_L + Kj_) + AsiS2q_i 






tmV2{x' - x) f ^ ei ■ T^ 5,,^, <5.,., + '^ 5^", ■ T^ 5. 

\xx (1 — a;)(l — x') 



S4S2 "S1S3 



8 (5siS3 ^: 



S3 ^S2S4, 



"S4S2 "S1S3 



, ni 



2m^ -^ -^ = Sqt . 

xx' {I - x){l - x') ^' 

The term No contributes to V in the order 



1 



m {iql - siql) i 1 J + q^ {ipl - Sip\) + -S2qz{q\ - isiql) 

m {iql - S2q\) ( 1 + ^^ ^ j - qz {ipl - Sspl) - -Siqziql - is2ql] 



,/\2 



2m' 



xx'{l — a;)(l — x') 



32m^ . 



(5.35) 



In these formulas we have used [/t_|_, k±] = e^jkIkI, Ey = Eijs and e* = — 4=(s,-i); also the 
following variables have been introduced 



q± = Kj_- K± , 

P± = ^ • 



(± = x,y) , qz=Pz- Pz 



(5.36) 
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We do not analyse here the expressions for A^i and A^2) where also in this form some terms can 
be identified as spin-orbit and spin-spin interactions in the transverse plane and in longitudinal 
(z) direction. 



Instead we follow |^, where an analogous calculation of singlet-triplet ground state mass 
splitting of positronium was performed in the similarity renormalization scheme. We therefore 
can drop in Ni, except for the leading order term 0(1), the diagonal part in spin space. Also 
the terms of the type / 
splitting, since 



x,y 
^_L Pz ) 



K_L Pz 5 n^Kj_ do not contribute to the ground state mass 



/ 



(iVV$Ioo(p) 4 "^100(^0 



(5.37) 



average over directions, gives zero. 

In the leading order 0(1) of nonrelativistic expansion we obtain the following ee-potential 

1 11/ 



V'('>\p',f) 



2(27r)3 4m2m 



X 



(p 



It 






IQe m 



m" 



ii 






^SlS-i ^S2Si 



a 



27r2g2 
V{t) f 



1 + 



1 



p2 

2^^ 

^\ 

2m?] 



"S1S3 "S2S4 



(5.38) 



where q = p' —p, and we have done in the last expression the Fourier transformation with 
respect to g*to the coordinate space. In the leading order of nonrelativistic expansion we have 
reproduced the Coulomb potential, defined before as the leading order term in BSPT. 

We combine this expression with the kinetic term from the Schrodinger equation, eq. (|5.13|) , 
and write it in the form 



1 L I V{r) 
m \ 2m 



p^ + V{r) . 



(5.39) 



Here the potential V{r) plays a different role in the two terms. In the first term, corresponding 
to kinetic energy, it generates an effective mass of the electron, which depends on the relative 
position and manifests the non-locality of the interaction. The second term is the usual potential 
energy, in our case, the Coulomb interaction. 

The energy of the Coulomb level for positronium with quantum numbers {nlm) is the 
standard one 



5W =< $„;^|y'(°)|$.,„ >= / d'p<fp'^l,^{p) 
eq. (|5.40|) the following representation 



a 



1 



^nlmip') 



ma 



, (5.40) 



2'K'^{p — p'Y) Ar? 

where the Coulomb wave function ^^nim is defined in eq. ( p.l4| ),eq. (|5.15|) . We have used in 



{P-P'f 
1 



[u — u 



'\2 



d^p = dVt^ 



E— >;(fip)>;*(fip' 

' el + f V 
2e„ / 



[u — u 



l\2 



(5.41) 



A7. 



and also orthogonality of the hyperspherical harmonics 

'dnY;Y^,=6,,,. (5.42) 

More details can be found in [^. 

Now the result of the first and the second order bound state perturbation theory for positro- 
nium ground state [n = 1) is presented. These are the corrections to the leading order binding 
energy B^' eq. ( |5.40| ), defined by the nonrelativistic Schrodinger equation for positronium 
eq. (HI). 

The next to leading order O (^) 



5V^'^ 



1 1 f e^' 



2(27r)3 4m2m V ^V 

X {8m{iql - Siql)6siS3^s2S4 - 8m{iql - S2qi)Ss^s3^s2S4) (5-43) 



contributes (because of the spin structure) to the second order of BSPT: 



S'^'B^ E <^^■°°l^^"'l^^^^^^><tr^'^^"''^''°°^ (5.44) 

-^-^ r(0) n(0) V J 



Recall, that /i = {n,l,m), the usual principal and angular momentum quantum numbers of 
nonrelativistic positronium. The order O i (^j 1 (cf. remark after eq. ( ^.36| )) is 

1 1 1 / g^ \ 

^^^^^ " 2(27)3 4^2^ {8(^^^^s^s2Ss,s■,Ss2S4 + 8e'^6s,s2SsssJs^ss + 4:e^Ss,s;Ss2s, j (5.45) 

and contributes to the first order of BSPT: 

d^B =< $100 1 (51^^^^ I $100 > • (5.46) 



Both contributions were calculated in [^ with the result 

5B = 5^'^B + 6^^^B 
5 

< 115511 > = ma^ 

' ' 12 

<2\SB\2> = <3\6B\3>=<4\SBU>=-ma\ (5.47) 

6 

where the eigenvectors in spin space are defined as follows: 

11 >=^(l+ ->-!-+», 

|2>= ^ (!+->+!-+>) , |3 >=!-->, |4 >=!++>. (5.48) 

V2 

Using the relation between Coulomb energy units and Ryd = ^ma"^, we get the standard result 
for the singlet-triplet mass splitting of positronium 

7 
Btripiet - B singlet = g^^Ryd + 0{ma^) (5.49) 

The degeneracy of the triplet ground state n = 1 signals the rotational invariance, that is 
not-manifest symmetry on the light-front. 

4.'^ 



Brisudova and Perry [^] tried to bring the effective light-front Hamiltonian of positronium, 
obtained in the second order in couphng by similarity transformation, to the rotational invariant 
form. They succeeded to get the correct spin-spin interactions, namely they obtained the 
familiar Breit-Fermi spin-spin and tensor terms. They failed to reproduce the standard (equal- 
time) spin-orbit interaction using the effective Hamiltonian in the order O(e^). The reason can 
be the following. 

Each spin enters the electron-positron interaction with a factor of order q/m as compared 
to the leading Coulomb interaction, that is of order g~^. Thus the two-spin interaction enters 
only in the order g° or higher. The only contribution to fine structure splitting, a^, comes from 
order e^ (from the terms e^g°). The same holds for the spin-triplet splitting (which is quadratic 
in the spin), discussed above. The contributions to the spin-orbit coupling are of the order q~^. 
In the order a^ also term e^g~^ is important, its contribution to the mass spectrum in order a^ 
must be considered together with the contribution of e^g~^ term. To get the correct spin-orbit 
splittings one has to derive the effective interaction to the order e^. 

Spin independent interaction is of order g~^ and contain also subleading terms in q (g~^, g°). 
The spin-independent term e^g° contribute to the fine structure splitting, in the order a^. In 
this order also terms of order e^q~^ and e^g~^ are important. Thus to obtain all contributions 
of order a^ one has to consider the effective interaction in order e^, e^ and e^. 
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Chapter 6 

Positronium spectrum (numerically) 



This part of the work is done in collaboration with Dr. U.Trittmann. In this chapter we suggest 
the form of the effective electron-positron interaction, that preserves the rotational symmetry 
(at least on the level of the mass spectrum). Using this interaction the light-front integral 
equation for positronium bound states is solved numerically^. 

6.1 Light-front bound state equation 

The effective light-front eigenvalue equation for positronium bound states reads (see subsection 

EH) 

H'ff\{ee)n>=M^\{ee)n>. (6.1) 



where n labels all quantum numbers, and the effective light-front Hamiltonian consists of the 
free (noninteracting) part and the effective electron-positron interaction 

H-ff = JjiO) + yeff (g 2) 

Note, that the rescaled value (that does not depend on P"*") of the effective electron-positron 
interaction, obtained in the chapter ^, stands in eq. (|6.1|) ,eq. (|6.2|) . The integral light-front 
equation in the momentum space, corresponding to eq. ( |6.1| ), is given 

= Zl / 0/0 \t < ^' ^-L' -^1' >^2\V''^^\x, k'j^; X[, A2 > ipnix', fcl; A'-^, A^) (6.3) 
i— ', Jd 2[27t)'^ 

where the wave function is normalized 

rjncd r\ i — * — * 

Y. 2f27r)3 ^"'-^^ ^^' ^^' ^2)V'n'(a;, ^±; Al, A2) = 5nn' (6.4) 

Ai,A2 ^ ' 

For practical purposes we have chosen Jacobi momenta as depicted on fig. (4). The integration 
domain D is restricted by the covariant cutoff condition [^5] 

^ +k\ ^ .2 , ._2 



x{l — x) 



< A^ + 4m^ (6.5) 



^ All numerical calculations were performed by Dr. U.Trittmann. The calculations are preliminary. 



which allows for states to have the kinetic energy below the bare cutoff A. 

In order to introduce the spectroscopic notation for positronium mass spectrum we integrate 
out the angular degree of freedom y?, where k± = /cj_(cos(/9,sin(/9), by substituting it with the 
discrete quantum number J^ = n, ra G Z (actually for the annihilation channel only | J2I < 1 is 
possible) 

< X, kr, Jz, Ai, A2|V"^^^|x', k'j_] J'^, X[, A2 > 

= ^ r^y^e-^^^ f'''d^'e^LW^_^^^)^^^k^^^.x,,X2\V^"\x',k'^,ip';X[,X', > 
In Jo Jo Z[Zn)'^ 

(6.6) 

where Lz = Jz — Sz] Sz = ^ + ^ and the states can be classified (strictly speaking only for 
rotationally invariant systems, that is the case in the nonrelativistic limit considered in the 
chapter |^) according to their quantum numbers of total angular momentum J, orbit angular 
momentum L, and total spin 5*. It should be noted that the definition of angular momen- 
tum operators in light-front dynamics is problematic because they include the interaction (see 
introduction in the chapter ^. 

We proceed now to solve for the positronium spectrum in all sectors of Jz- We formulate 
therefore the light-front integral equation eq. (|6.3| ) in the form where the integral kernel is given 
by the effective interaction for the total momentum Jz eq. (|6.6| ). After the change of variables 
eq. ( |5l^) {k±;x) = {k±,ip;x) — > p = {k±,kz) = (/i sin 6* cos ip, ^sinOsmip, ^ cos 9) 

where the Jacobian reads 

/i^ m^ + /x^(l-cos^g) 

'^ = ^ — r^^ — T^2 — si^^ ^-s 

one has the following integral equation 

+1 ,,/i'2m2 + /i'2(l-cos2^') 



{M^-4:{m'^ + lx^))i!nil-i,cose;Jz,XiA2)+ Y. [dft'l dcos9' 



2 (m2 + /i'2)3/2 

X < /i, COS 6*; J^, Ai, X2\V^^^\iJi! , cos6'; J'^, X[, X'2 > 'ipnifJ'', cos 6*'; J^, A'^, Aj) = 



(6.9) 



The integration domain D eq. (|6.5| ) is given now by ;U G [0; ^]. Neither Lz nor Sz are good 
quantum numbers; therefore we set Lz = Jz — Sz- 
The wave function is normalized 

J2 djj,dcos6ijj^{iJ,,cos6;Jz,Xi,X2)'ijjn'{fJ',cos6;Jz,Xi,X2) = 6nn' (6.10) 

Jz, Xl,X2 

where n labels all quantum numbers. 



The integral equation eq. ( |6.9| ) is used further to calculate the positronium mass spectrum 
numerically. 
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6.2 Brodsky-Lepage prescription of the light-front dy- 
namics: effective electron-positron interaction 

In this section we write the effective electron-positron interaction, using the hght-front prescrip- 
tion as formulated by Brodsky and Lepage [E3]. The problem in the light-front field theory 
is, that parity and rotational invariance are not manifest symmetries, leaving the possibility 
that approximations or incorrect renormalization (in our case not gauge invariant regulariza- 
tion) might lead to violations of these symmetries for physical observables. The formulation 
for light-front field theory of Brodsky, Lepage is for practical use, since it enables to pick up 
the rotation violating part in the electron-positron interaction and to cancel its contribution in 
the mass spectrum. 

In what follows we write the second order solution of the flow equations in a useful for 
our calculations form (subsection 6.2.1) and exploit it below to calculate the effective electron- 
positron interaction, using Brodsky-Lepage prescription of light-front QED (subsection 6.2.2). 

6.2.1 First and second order solutions of the flow equations 

In this section we express the solutions of the flow equations in the first and second orders of 
perturbation theory through the similarity function. The similarity function defines the 'rate' 
how fast the 'rest' sectors are eliminated with the flow parameter / (or effective UV-cutoff A); 
exact definition is given below. 

We consider the flow equations, written in the form eq. (|2.31 ) 



TT = [^7' ^d + Hr\ij — {Ei — Ej)[Hd, Hr]ij H 

dl dl Uij 

,,^.i^,.ff.,], + _l_(_^^) (6,1) 

where the following conditions on the cutoff function in 'diagonal' and 'rest' sectors, resp., are 
imposed 

Udij = 1 

Urij = Uij (6.12) 

In the perturbative frame we break the Hamiltonian as 

i7 = i7od + E(^i"^ + ^r^"^) (6-13) 

n 

where H^""^ ~ e", e is the bare coupling constant (here we do not refer to the definite field 
theory). To the order of n in coupling constant the flow equations in both sectors are given 

r(") 



dH. 



''^ -Y.W^\Hi''-%^, 



dl 

-^ = Uv'''. Hf' + Hi-%., - (E, - E,) E[^f , Hi--%,, ^d_^H^ 

k k '^3 

M _yr^,rT{k) TT{n-k)] , ^ | ^'^ij ^rij 

"^"^ ~Y' '^ ^"'^ E,-E, [ dl U,, 
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We solve these equations in the leading order for the Fock state conserving sector. 



dH. 



dl 

Irij 



(2) 



(1) _ 1 dH^^j 

Ei — Ej dl 



dHl^_,lu„Hl^ 



dl dl Uij 



(6.15) 



and if^jj = -f^^ji- Explicitly one has 

dl ^[E.-Ek dl ^^'^ Ej - Ek dl 



''"^^^^~-^"^^^~^^ fJl = Q) 



Hllil)=K]{l = 0)-^^i^ (6.16) 



where we have introduced the function fij defining the leading order solution for the 'rest' part. 
Further we refer to it as similarity function. Here 

/,,■(/) = «,,(/) =e-(^'-^^-)'' (6.17) 

The similarity function /a (A) has the same behavior (when A ^ oo /a (A) = 1, and when 
A — > /a (A) = 0) as the cutoff function ma(A). 
Making use of the connection / = 1/A^, we get 



dX Y V^* ~ ^'^ dX "^^ Ej - Ek dX 

HiiliX) = H^iA -. oo)— M^ (6.18) 

Neglecting the dependence of the energy Ei on the cutoff, one has 

^g(A) = i^g(A - oo) + 5: {H^liA - oo)/7« (A - oo))^ 

where A is the bare cutoff; the sum J2k is over all intermediate states; and the label 'd' denotes 
the 'diagonal' sector. 

In the case of other unitary transformations, eq. ( |2.33|) ,eq. (|2.34|) , the similarity functions 
/(A) are given 

— -^ = Uij[r], Hd + Hr\ij + Tij—^—^ 
dX dX u 



V 



n., = -^UH, + H: ^""-^^ 



E.^\^^^^^^^^^-^^, 



/,,(A) = «.,(A)e^-W 

Ui, + n, = 1 (6.20) 
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and 



"^^'^ -u.Av,H, + Hr]., + '^''''^'' 



d\ 



dX u 



V 



Vij 



Ei — Ej 



rij[r],Hd + Hr]ij - '^ '^ 



d\ u 



^3 



fij{\) = Uij{\) 



(6.21^ 



where u{X) is the cutoff function. One can choose Uij = 9{\ — |Ajj|), where Aj^ = I]fcii -Si.fc — 

fc=l ^j,k- 

The equations eq. ( |6.19| ) are the same for all unitary transformations, given above up to the 
choice of the similarity function /(A). Specifying the function /(A) we get the interaction if^ , 
generated by different unitary transformations. 



6.2.2 Effective electron-positron interaction 

The exchange channel brings the dominant contribution to the mass spectrum. We focus 
therefore on the electron-positron interaction in the exchange channel. We use eq. ( |6.19| ) to 
calculate the generated electron-positron interaction in the second order in e. In this case 
the 'diagonal' sector (denoted 'd') is |ee > sector, and the matrix element is given fig. (4) 
< pi, Ai;p2, A2I...IP1, Ai;p2) A2 >; H^^\K -^ 00) is the electron-photon coupling term with the 
bare coupling constant e; the initial value of generated interaction is given iJ^ (A -^ 00) = 0. 
The generated interaction is given 



ygen ^ _^2 ^ ^^.^u 



ex{D,,D2) = J^ 




(6.22) 



where we sum {J2k iii eq. (|6.19| )) the two terms corresponding to the two time-ordered diagrams 
with g"*" > and g"*" < 0; A is the 'running' cutoff, that defines the continuous step of the unitary 
transformation; /a (A) is the similarity function, arising from the unitary transformation; the 
function 9e{q^) restricts the longitudinal momentum of intermediate photon. The similarity 
function /a(A) and the cutoff function Os^q^) are specified below.0 D^^i^q) = J2x ^^(A, g)e*(A, q) 
is the polarization sum; the null vector i]^ = (0,r7+,0,0), ri^rj^ = is given below; the energy 
denominators in the exchange channel are given (fig. (4)) Di = pf — Pi — q~ and D2 = 

P2 ~ P2 ~ Q] Q = Pi ~ Pi is the exchanged photon momentum, with q~ = ^. The notation 
< 'j'^'j'^ > is introduced for the current- current term. In the exchange channel this matrix 
element is given fig. (4) 



< TT > \exch 



u{piAi)_ ,, n(Pi,Ai)i;(p2,A2) ^ v{p2A2) p+2 



(6.23) 



^The similarity function /a plays the role of transverse UV regulator (see the chapter 0), and the function 6,. 
regulates the longitudinal IR divergences. On a tree level there is no UV divergences in four fermion interaction, 
but the longitudinal IR divergences are present. 



4P 



where piiP^ are light-front three-momenta carried by the constituents, Aj, \ are their hght-front 
hehcities, u(pi, Ai), f (p2, A2) are their spinors (see below); index i = 1,2 refers to electron and 
positron, respectively; P = {P~^,P-^) is light-front positronium momentum. 

Below we use the light-front conventions formulated by Lepage and Brodsky EH] (see also 



1 2^). The polarization sum is given 



(6.24) 



where the metric tensor(s) 
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and the null vector is rj'^ = (0, 2, 0). The Dirac spinors are given 

X(T), forA = +l, 



\ 




-1/ 



m(p. A) 
v{p, A) 



1 

7p 

1 



/3m + d±'p± ) X 



Pm + a±_p±_ ) X 



X(i), forA = -l, 
Jxd), forA = +l, 



/p-r y / lx(T), for A 

Here /? = 7*^, a = 7*^7; and the two x-spinors are 



x(T) 



1 



V2 





1 
VO/ 



and 



X(i) 



v^ 



1 


V-1/ 



(6.25) 



(6.26) 
(6.27) 



(6.28) 



These conventions are used in Appendix ^ to calculate the matrix elements of the effective 
interactions 0. 

Using the symmetry 



fx{-D) = h{D) 



(6.31) 



■^ In the chaptersplandH we have used the prescription of the hght-front field theory as formulated by Zhang 
and Harindranath p4| . They have used the following conventions: the polarization sum is 



d^^Aq) 



q 



1 



+2 '^M^i' + —inf^Qu + fli^Qf, ) - 9f,u 



(6.29) 



with the null vector 77^ = (0,1,0); the four-component spinors w(p. A), v{p,X) are given through the two- 
component spinors xa in eq. (??) and eq. ( 3.10 ), chapter]^. 
Using the above equations, we get 



<7^7">|e.ch^M-(9) = M, 



ex 
2ii 



(6.30) 



Jex) 



where M^^j' = [x^.T'i^p'^^pi^p'^ -pi)xA'J [x^' r^(-P2, -P2> -(p'l -Pi))Xa2] i^ the matrix element, obtained in 



the chapter ^ eq. ( |4.29| ) direct from the two-component field theory eq. (3.3)-eq. (3.7). We reproduce then the 
form of generated interaction, obtained in the chapter ||eq. (4.3C). 
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we have the following generated interaction in the electron-positron sector 



(6.32) 
We specify the cutoff function 6'e in a form [^| 

e, = eiq+-e)F{q+,q^) (6.33) 

2 

where e is 'small', say £ ^ ^ and A is the bare cutoff. For our purposes we do not need the 
explicit form of the function F{q~^, qj_), it defines the upper boundary for Ig+I. The combination 
{Oe{q~^) + (^e{~q~^)) restricts |g^| to be above e. In the integral in dq'^ this ensures the symmetric 
cutoff for longitudinal photon momentum 

/— £ r+oo 

dq+ + / dq+ (6.34) 

00 Je 

The resulting generated interaction is given at the value A = 0, V^^e" _ V^^q. To get the 
effective electron-positron interaction in the light-front dynamics we sum the resulting generated 
interaction and instantaneous photon exchange arising in the light-front QED 

yeff ^ ygen ^ yinst ^g 35^, 

where each term is defined 



^3- = -e' < YY > -D,M) ' " + 



1 

2' 



e(A,/^2) = Qx=,{D,,D2) = r ^^^fx'iD2)dX' (6.37) 

JO dA' 

For both interactions the prescription eq. ( 6.33|) (eq. ( |6.34|) ) to treat infrared divergences, small 
q~^, is imposed. We combine all the terms together with the result 

This equation does not depend on the explicit form of 0-factor. The generated interaction 
has two types of infrared singularities: -^ and ^ types. We see that the -^ singularity of 
generated term is cancelled exactly by the instantaneous term in the effective electron-positron 
interaction. Further we show that the cutoff condition eq. (|6.33|) ensures the cancellation of 



-^ type singularity in the physical observables, calculated from the effective electron-positron 
interaction. 

We introduce 

2d = Di-D2; D = £l±£l (6.39) 

that gives 

Di = D + d; D2 = D-d (6.40) 



As g^ — i> we expand the effective interaction eq. (|6.38|) in the series in terms of -^ << 1 
{d^{A^Q-Ml)/P+-D^{-q~)) 



q+D 

1 

+2 



1 - - [0(Di, D2) - 6(^2, D^)] + (— ) - (— ) [0(Di, D,) - 6(^2, A)] 
1 



q^ 
X (-;^[epi, D,) - 0(^2, 1^1)] + i^,) - (|^)[0Pi, 1^2) - e{D„ Di)]) 

= n^'' + EAV;«+EAna. (6-41) 

where we have used the identities 



Di D2 2 V/^i £'2/ 2 VA £'2/ 

e(Di,D2) + e(D2,A) = i (6.42) 

index [i) denotes the order with respect of -^. The leading order term is given 

Vo" = e'< YY > 9,.4t^ (6.43) 

q+U 

In the leading order of nonrelativistic approximation \p\/m << 1 this term eq. (|6.43|) gives the 



3-dimensional Coulomb interaction [53| 



K" - -'-^ (6.44) 

q^ 

where q(qz,(l±) = p' — p is the exchanged momentum. Hence the effective electron-positron 
interaction eq. ( |6.38| ) produces Bohr energy levels [Q. The nonrelativistic expansion of the 
term eq. ( |6.43| ) up to the second order O (^j gives the familiar Breit-Fermi spin-spin and 
tensor interactions |^, that insures the correct spin-splittings for the positronium ground 



state []39 |. 



Corrections AV^*^^ ^^d AVj|*^ arise due to the unitary transformation performed, i.e. that 
are the corrections due to the energy denominators in the ^^g'^^ term and the "'/^^r]" term. The 



first order corrections 0{d/D) were estimated in [^,||55|, using the explicit form for similarity 
function fx{D). We are interested here in general properties (independent on the choice of 
/a(-D)) of the effective interaction, particularly in the origin of infrared divergences. 

In eq. ( |6.41| ) the product oi"r]^r]'l term with the term ~ (-§) gives the singularity of -\ type 
(when {Q{Di, D2) — Q{D2,Di)) ~ const with respect to (d/D)). To find the corresponding 
infrared counterterm we integrate the singular term over 'external legs' (that correspond in 
QED to physical particles) with the wave packet function, J dq'^d'^q±(f){p + q), where p ~ Pi,P2- 
The integration of -^ type term leads to logarithmic infrared divergences, that are cancelled 
from small positive and negative longitudinal momenta due to the symmetric cutoff condition 
eq. (|6.34|) . Therefore, there are no infrared counterterms to be introduced in the second order 



O(e^) in the effective interaction. Moreover, this is true in the order O(e^) for all tree level 
diagrams in QED. 

The same cancellation of the infrared divergent contribution from "rj^ri'l term occurs, due 
to the symmetric cutoff, in the spectrum of masses for positronium. 

What about the finite corrections from "rj^rjl term? The leading order finite correction 
from "rf^ri'l^ term eq. ( 6.41|) is ~ f]^Vp-T2{j)Y^ which is of the order e^g". The interaction of 



light-front field theory V^'^ eq. ( |6.52| ) and the effective interaction V^^^ eq. ( |6.38| ), generated 



by the unitary transformation, have both a leading Coulomb behavior, but they differ by spin- 
independent "Tj^rj'l term (excluding the divergent part) in the order e^g°, which contributes 
to the mass in the order a^. In the order of fine structure splitting a^ also terms of order 
e^q~^ and e^g~^ are important. These terms arise from the next orders transformation of both 
electron-photon vertex and also instantaneous term. Since the Coulomb interaction, that is 
the only spin-independent part of electron-positron interaction, arise from " g'' term in the 
effective interaction V^^^ , we expect that the spin-independent "rj^ri'l term in the order e^g° 
will be compensated in the mass spectrum by the corresponding terms of the order e^ and e^. 
Compare effective electron-positron interaction eq. ( p.38| ), generated in the second order in 
e by flow equations, with the electron-positron interaction arising from perturbative photon 
exchange. The electron-positron interaction in the light-front perturbation theory is given by 
a sum of dynamical photon exchange V^p'*"* and the instantaneous interaction, V^^ = V^^"^ + 
yznst^ The instantaneous term is given by eq. (|6.37 ), and the V^^"^ term, defined by the diagram 



with one photon exchange, is |45 



VP'^ot ^ _^2 < ^,^. ^ f Milz)^^^(^) 1 + ^Ajllo^^i-q)^] (6.45) 




where the energy denominators are given D± = J2incP~ ~ J2mtermP~ ^ with the sums over the 
light-front energies, p~ , of incident (inc) and intermediate (interm) particles; and ± denotes 
two different time orderings of the photon. For the process depicted on fig. (4) we have D^ = 
P^ — Pi — P2 — (T and D_ = P — pi — P2 + Q^ , where q = p'l — pi, P = M^/P^ (let 
P^ = 0). One can rewrite eq. ( |6.45|) in the form [^ (also see Appendix |^) 

yp'^ot ^ _^2 ^ ^,^u > Msll±MzSllD^^(q)i (6.46) 



where 



9,iq+)D+ = -^,(-g+)D_ = D (6.47) 



Combining instantaneous eq. (|6.37|) and one photon exchange eq. ( |6.46| ) terms together, we 
have 

V^^ = e'< YY > 9,u^ -e'< YY > V,Vu^ U " ^) (6-48) 

q+D q^-^ \ DJ 



where the prescription eq. ( |6.33| ) (eq. ( |6.34] )) to treat infrared divergences is imposed. 



The "rj^rj'l^ term carries the infrared divergent part of the interaction V^'^ . We approximate 
the positronium mass as average over the masses of initial and final free (not bound) states 

Ml = -^— — ^ (6.49) 



F,?, 



where pr + P2 = M^/P+, p'f + p^ = M'q^/P+, P = M^/P+, and P{P+,Pi_ = 0) is the 
positronium momentum. In this approximation the energy denominator of perturbative photon 
exchange diagram is 

D^D = ^^±^ (6.50) 

It is obvious, when one notes (see also Appendix |B]) 

D1 + D2 p'f -Pi +P2 - P2 - Pi + P2 +Pi +P2 

D = — ^ — = ^ q = ^ Pi -P2 -Q 

+ P2 +Pi +P2 '- - , -\ (a KA\ 
n Pi -P2+Q (6-51) 



Given the eq. (|6.5CI|) , the instantaneous interaction \/*"^* eq. (|6.37|) , which is the source of 
infrared divergences, is precisely canceled by the part in V^^°^ corresponding to emission and 
absorption of longitudinal photons. The resulting electron-positron interaction, obtained in the 
light-front perturbation theory, is given in this approximation eq. (|6.49|) 

V^ = e'< 7^^ > g,.^ (6.52) 



By the special choice of similarity function this answer eq. ( |6.52|) can be generated by the flow 
equations. When the similarity function is 



DP' 
\ 



fxiD) - 

D' = signD (6.53) 

the 6-factor eq. ( |6.37|) is 61 = -jj-t^, and "rj^rj'l term is equal to zero in eq. ( |6.38|) . The 



effective electron-positron interaction eq. ( |6.38| ) is given by the result of perturbation theory 
eq. (PI) 

y'"\MD) = V^ (6.54) 



where the condition eq. ( |6.49| ) (eq. ( |6.5CI| )) is imposed. This is the remarkable result, that shows 



that specifying the unitary transformation in a proper way one can get rid of "rjf^rj'l term in the 
effective interaction eq. (|6.38|) , the term that causes the infrared divergences in longitudinal 
direction. In all other cases we need to introduce the symmetric cutoff condition eq. ( 6.33| ) 
(eq. ( |6.34| )) on g"*" to cancel the divergent contribution in physical observables. 

What is the status of "77^77" term in the effective electron-positron interaction V^^^ eq. (|6.38|) ? 
The "rjfj^'q'l term describes instant emission and absorption of 'longitudinal' photon and is specific 
in the light-front gauge computations. As was discussed above, this term is not desirable in the 
electron-positron interaction and can be considered as a consequence of the unitary transfor- 
mation performed. The physical reason for its appearance is the violation of Lorenz and gauge 
symmetries by the derivation of the effective, renormalized Hamiltonian. The first problem in 
the light-front field theory is that whenever the generator of a symmetry is dynamical (contains 
interactions) it is somewhere between very difficult and impossible to monitor and maintain 
that symmetry at each step of a calculation - unless of course one can solve the theory ex- 
actly. This concerns parity and rotational invariance, that are not manifest symmetries on the 
light-front. The second problem is that whenever we use the Hamiltonian technique, the naive 
regularization by introducing the cutoffs breaks the gauge invariance (and also Lorenz covari- 
ance), and forces the bare Hamiltonian to contain a larger than normal suite of counterterms 

H4 



to enable a finite limit as the cutoffs are removed. One way is, that the counterterms are then 
adjusted to reproduce physical observables and to restore the symmetries broken by cutoffs. 
The other way is to find the gauge invariant procedure for regularization of Hamiltonians. The 



attempt in the latter direction was made by Brodsky, Hiller, and McCartor |^. Solving the 
Yukawa theory they have tried to preserve more symmetries by using Pauli-Villars procedure 
for regularization. 

In the next section we use the effective electron-positron interaction V^'^'^ eq. ( |6.38| ), namely 
its "g" part, to calculate mass spectrum and wave functions for positronium. 

6.3 Mass spectrum and wave functions of positronium 

We start with the effective electron-positron interaction, generated in the second order in cou- 
pling by flow equations, eq. ( |6.38| ) 

- e^ < YY > ,l^,(D. - D..) {^^ - 5^) (6.55) 



where we impose the prescription eq. ( |6.33| ) (eq. ( |6.34| )) to treat infrared divergences, small q^ . 
We use Jacobi momenta, depicted on fig. (4) 

Pi(a:P+,xPi + fc±) 

P2{{l-x)P+,{l-x)P^-k^) (6.56) 

and corresponding for the momenta Pi,P2) ^^^^^ ^ i^ the light-front fraction of electron momen- 
tum and P{P^, P±) is the total momentum of positronium. For convenience we introduce 

^1 = 7^ = -^' ^^=PT = -^ (6-57) 

and from now on we use the rescaled value of the cutoff A —>■ X'^/P^. The effective electron- 
positron interaction is written then 



V Ai A 



2 



e' < ^^Y > Vt^VuTTT^i^i - ^ 



1 ,. ~ ,/e(Ai,A2) e(A2,Ai: 



2g+2^ ^ ^^V Ai A2 ; 

e(Ai,A2) = /f ^^^/A'(^2)rfA'2 (6.58) 

where the energy denominators in eq. ( |6.58|) read 



- {xk'j_ — x'k±y + m'^{x — x')'^ - ~ 

Ai = ; A2 = Ai|a;.^(i_a;)^^/-^(l_a;/) 

Ai = -^ ; A2 = -^ (6.59) 

X ~ X X — X 

Note Ai, A2 are positive definite. 



To perform the angular integration eq. (|6.6| ) analytically we choose the similarity function 
as follows 

/a(A)=m,(A) = 0(A2-|A|) (6.60) 

where for simplicity we use the sharp cutoff function u\{A). Then the effective interaction 
reads 

where we have used 9{x) — 9{—x) = sign{x) 0. 

The general matrix elements for the effective interaction eq. (|6.61|) depending on the angles 
(p, ip' (actually on the difference {ip — ip')) < x, k±, ip; Ai, X2\V'^'^^\x', k'j_, ip'; X[, X'2 > and also the 
matrix elements of the effective interaction after the angular integration eq. (|6.6| ) for the total 
momentum J^ < x, k±; J^, Xi, X2\V'^-^-^\x', k'j_; J^, A'^, A2 > are given in Appendix ^. 

The positronium spectrum is calculated numerically, using the integral equation eq. ( |6.9| ) 
with the matrix elements given in Appendix |B[ We use for the numerical integration the 

"* In this particular case, when the similarity function is given by eq. ( 6.6C ), the divergence arising is of the 
vents from the first glance the cancellation of logarithmic divergences in the integral 

dq+ (6.62) 



9+l>e 

We have chosen the sharp cutoff for simplicity, but generally the smooth cutoffs are preferred to avoid non- 
analyticities in the structure of counterterms pO| . We demand only, that the cutoff function u{x) with x = 
■LJ- — > W- is 1 for small arguments and vanishes for large arguments. Explicit one can choose 

u(x) = I for < X < — 
w J - - 3 

1 2 

u{x) monoton falls from 1 to for — < x < — 

2 
u{x) =0 for - <x<l (6.63) 

o 



and to simplify the calculation of the integral eq. (6.58), we assume 7"^ ~ (5(1 — x). Then, when the similarity 
function (the function that defines the first order solution in the perturbation theory of the flow equation) is 
equal to the cutoff function f{x) = u(x), the effective electron-positron interaction reads 



In this case the prescription eq. (6.62) ensures the cancellation of logarithmic divergences, and one does not need 
to introduce counterterms. In the integral equation eq. (6.3) we have the following restriction of the integration 
domain 

'' dx' ^ I d{x'-x) (6.65) 



resulting, that the divergent part in "rj^rj'l term of the effective interaction eq. (6.64) does not contribute to 
mass spectrum. 



Gauss- Legendre algorithm (Gaussian quadratures). To improve the numerical convergence the 
technique of Coulomb counterterms is included. The problem has been solved for all components 
of the total angular momentum, J^. Since we calculate the values of an invariant mass squared, 
a strong coupling constant a = 0, 3 has been chosen. The latter means also, that the method 
of flow equations is applicable in non-perturbative regime. We get the ionization threshold 
at M^ ~ 4m^, the Bohr spectrum, and what is more important, the fine structure. The 
agreement is quantitative (for the lowest eigenvalues), particularly for the physical value of the 
fine structure constant a = j^. 

The figure (5) is the summary of the spectra for different components of the total angular 
momentum, J^. As one can see, certain mass eigenvalues at J^ = are degenerate with certain 
eigenvalues at other J^ to a very high degree of numerical precision. As an example, consider 
the second lowest eigenvalue for J^ = 0. It is degenerate with the lowest eigenvalue for J^ = ±1, 
and can thus be classified as a member of the triplet with J = 1. Correspondingly, the lowest 
eigenvalue for J^ = having no companion can be classified as the singlet state with J = 0. 
Quite in general one can interpret 2 J^ ^ax + 1 degenerate multiplets as members of a state with 
total angular momentum J. One can get the quantum number of total angular momentum 
J from the number of degenerate states for a fixed eigenvalue M^. Such form of spectrum is 
driven by rotational invariance. It is a remarkable result, we restore the rotational symmetry 
in the light-front positronium calculations on the level of mass spectrum. 

The integral equation is approximated by Gaussian quadratures, and the results are studied 
as a function of the number of integration points A^, as displayed in Figures (6) and (7). On 
the fig. (6) the mass squared eigenvalues M^ for J^ = are shown as functions of the number 
of integration points N = Ni = N2. One sees, that the results stabilize themselves quickly. 
To find out if the degeneracy obtained is merely a numerical artifact, or a property of the 
positronium model, consider fig. (7). The mass (squared) discrepancy between the J^ = and 
Jz = 1 eigenvalues is plotted versus the number of integration points N for three different 
states. It is important, that the convergence of the value AM^(l^S'i) with N occurs. But the 
mass gap does not converge to zero, rather to some small value ~ 10~^. The mass gap for 
2^Pi state vanishes as A^ increases, while for 2^Pi state AM^ converges again to a finite small 
number ~ 10~^. Kaluza and Pirner |^ and Trittmann and Pauli [Q found that in light-front 



perturbation theory there is a discrepancy between the case of J^ = and J^ = 1. 



F,7 



Chapter 7 

Renormalization in the light-front QED 



7.1 Introduction. 

In this section we consider the renormahzation group problem for hght-front QED, arising first 
in the second order in couphng. Renormahzation in the hght-front field theory is much more 
complicated than in covariant Feynman perturbation theory, and is nontrivial already in the 
leading order for QED. Generally, there are three renormalization problems that are signaled 
by the divergence of the free energy: 

(1) Pi; 

(2) p+^O; 

(3) J9+ = 0. 

The divergences from (1) are 'ultraviolet' (UV), producing counterterms that are local 
in the transverse direction. While the counterterm structure is richer than what is required 
in manifestly covariant perturbation theory, in perturbative regime one has that masses are 
relevant operators, and other gauge interactions are marginal. 

The divergences from (2) are 'infrared' (IR), producing counterterms that are nonlocal at 
least in the longitudinal direction and which may also be nonlocal in the transverse direction. 
These infrared divergences do not appear in Feynman perturbation theory in covariant gauges. 
They must cancel perturbatively for gauge invariant matrix elements in QED. Coupling co- 
herence fixes the infrared divergent part of the Hamiltonian as a power series in the gauge 
coupling, allowing us to next consider re-summations that preserve cancellation of all infrared 
divergences. 

Finally, (3) is the vacuum problem, that we do not consider here. The light-front vacuum is 
'trivial', except for zero modes. The light-front gauge A~^ = has a singularity at longitudinal 
momentum p~^ = 0, i.e. the zero-modes p"*" = should be treated specially, that is called the 
problem of zero-modes. One can connect (3) and (2), but we ignore the problem of zero modes 
in QED (at least here). 

In light-front renormalization group transverse and longitudinal directions scale separately, 
and transverse scaling runs the cutoff. We assume transverse locality and identify relevant 
and marginal operators using power counting, which is valid in the perturbative regime. The 
longitudinal scale is not important for the perturbative classification of operators. This means, 
that entire functions of longitudinal momenta (their dimensionless ratios), including any longi- 
tudinal momentum scale introduced by cutoff, appear in the relevant and marginal operators. 
The simplest way to adjust the new operators is to fix broken by the cutoff covariance (Lorenz 
covariance) and gauge invariance in physical observables. Generally, infinite number of coun- 



teterms must be added to the Hamiltonian to restore these symmetries in observables, and to 
obtain finite results. As a result, the renormalizability and as a consequence predictive power 
of the theory seem to be problematic. 



Coupling coherence, developed by Perry and Wilson [|J, says how to deal in practice with 
the functions that appear in marginal and relevant light-front operators and how to reduce the 
number (infinite number) of couplings, that arise by renormalization scaling when the sym- 
metries of the theory (rotational and gauge invariance) are broken by cutoffs. We demand 
that these additional counterterms required to restore the symmetries run coherently with the 
"canonical renormalizable couplings" - that is we do not allow them to explicitly depend on the 
cutoff but only implicitly through the "canonical renormalizable couplings" dependences on the 
cutoff. In other words, coupling coherence provides the condition under which a finite num- 
ber of running variables determines the renormalization group trajectory of the renormalized 
Hamiltonian. Also it is conjectured that coupling coherence restores (or reveal) symmetries 
broken by the regulator or vacuum. 

Coupling coherent Hamiltonian written in terms of dimensionless couplings for A << A 
satisfies 

Hx = f/(A, A)Ha{A, ca, rriA, wa, c(eA, mA))f/+(A, A) -^ i/A(A, ex, rux, w{ex, mx), c{ex, rux)) 

(7.1) 

where U{X, A) is the unitary transformation defined by fiow equations. The additional require- 
ment to eq. ( rnp is that all dependent couplings (represented by 'c' in the argument of the 
Hamiltonians) vanish when the independent marginal couplings are set to zero. In eq. ( [7. 1|) , ca 
and rriA are independent dimensionless marginal and relevant couplings, respectively; wa rep- 
resents the infinite set of independent dimensionless irrelevant couplings; c{ex, mx) represents 
the infinite set of dependent dimensionless relevant, marginal and irrelevant couplings. More 
about classification of couplings see |H3|. 



The initial bare Hamiltonian i^A does not satisfy eq. ( [7.1|) a priori, its form changes under 
the action of transformation U{X, A).Ha must be adjusted until its form does not change. 
This "adjustment" is the process of renormalization. Coupling coherence is a highly nontrivial 
constraint on the theory and to date has only been solved perturbatively. Further we present 
the solution of eq. (|7. 1| ) for QED in the second order in coupling (we calculate the electron and 
photon 'mass' operators), that turns out to be simple because coupling constant does not run 
until the third order. 

In the second order in coupling the program we proceed to find renormalized electron and 
photon masses is the following: 

(1) regularize the canonical light-front QED Hamiltonian. In this case, we introduce the bare 
cutoff A, that regulates the divergences in transverse direction, and introduce the second order 
mass counterterm 5M^ . We start with the bare cutoff mass m\ = m? + 5M^ , where m is 
the mass term in free Hamiltonian Hq; 

(2) perform the unitary transformation f/(A, A). The second-order change in the Hamiltonian 
is given in eq. ( |4.14] ) (section [4.2] ), which gives in one-body sector the shift to the electron 



(photon) mass 

~ / < [r]^'\Hee^] > one-body dX' ~ -(5Sa " ^^Sa) (7.2) 

where r]^^' is the first order generator of transformation, defined by fiow equations as r]^^' = 
[Hq, Hee-y], Hee-y is the electrou-photou vertex operator, and 6T,x is proportional to the value of 
the integral in eq. ( [7.2|) at the point A and defines mass correction, namely self energy term. 



Physically, in the electron sector the energy scales A down to A from photon emission have been 
'integrated out' and placed in effective interactions (the corresponding in the photon sector); 
(3) renormalize the Hamiltonian, i.e. adjust the Hamiltonian that its form does not change 
under the unitary transformation eq. (|7.1| ) using coupling coherence. As a result we get the 
renormalized (in the given order) Hamiltonian H\. In this case, since the electron-photon 
coupling does not run until third order, the constraint eq. (|7.1| ) is simple to solve: to second 
order the self-energy must exactly reproduce itself with A ^ A 

[m^ - (5Sa - <5Sa) + 5Mf^ + 0(e^)] = [m" + SU'i'^ + 0(e^)]A^A 

= [m'-5T.x + 0{e^)] (7.3) 

this fixes the counterterm 

5Mf> = -5Sa + O(e^) (7.4) 

i.e. in the second order electron (photon) mass squared runs coherently with the cutoff according 



to eq. (|7. 1| ) as 

ml = m^ - (5Sa (7.5) 

(4) using the renormalized Hamiltonian H\ calculate physical observables, that due to the 
constraint of coupling coherence are cutoff and renormalization scale independent (are finite as 
A — »• oo and £ — > 0, where A is UV regulator in transverse direction and e is IR regulator in 
longitudinal direction (generally, it can be a set of regulators), and are independent of A or any 
other scale introduced by renormalization), and manifest the symmetries of the theory broken 
by cutoffs. We calculate the physical mass of electron (photon), that includes the running 
mass m\ from Hx and in the second order the perturbation theory corrections. In the electron 
sector these corrections arise from the low-energy photon emission, since H\ still has the photon 
emission interaction below A of the form f\ J (Px^dx~Hee^. (One has in the photon sector the 
corresponding interaction term, where the energy of ee pair is restricted by A.) It turns out, 
that the electron (photon) perturbative mass-squared shift is exactly equal to the self-energy 
term obtained in (3) 

6m^ = <5Sa + O(e^) (7.6) 

therefore the physical mass-squared is given 

^Ihys = ""^A + ^^"^ = m^ — (5Sa + (5Sa = rn^ (7.7) 

where m"^ is the renormalized electron mass-squared in the free Hamiltonian Hq. This result 
is to be expected, that shows that the renormalization procedure using flow equations and 
coupling coherence is successful (at least in the second order). 

All details of calculation can be found in two sections that follow. The renormalization 
scheme presented seems to work trivially simple in the second order. But there are difficulties, 
that one encounters in practical calculations. The integral in eq. ( [7.2| ) is infrared singular 
(in both electron and photon sectors), resulting in infrared divergent self-energy and hence 
infrared divergent running mass m\, with new types of divergences (the electron mass has 
both linear and logarithmic infrared divergences). This is in accordance with the discussion 
presented above. The longitudinal cutoff violates boost invariance and the mass operator, that 
is a function of a longitudinal momentum scale introduced by the cutoff, is required to restore 
this symmetry in physical results. Indeed, as shown by Perry |^, the divergent self-energy is 



fin 



exactly canceled by perturbative mixing with small-x photons (see eq. ( [7.6|) ,eq. (fTT])), resulting 
in finite mass. 

We try to preserve boost invariance, that is manifest symmetry on the light-front, in the 
second order QED calculations on the level of renormalized Hamiltonian. We take into account 
the diagrams arising from the normal ordering of instantaneous interactions. The bare light- 
front QED Hamiltonian contains both the instantaneous interactions and the electron-photon 
vertex, therefore instantaneous terms must accompany the latter also by scaling from the ener- 
gies A down to A eq. ( |7.1| ). We get the running electron and photon masses free from infrared 
divergences, moreover both electron and photon mass conterterms contain the known from the 
covariant perturbation theory types of divergences. Unfortunately, the electron wave function 
renormalization constant Z^ contains mixing UV and IR logarithmic divergences and also pure 
IR logarithmic divergence; the renormalization constant Z2 in photon sector is IR finite. It is 



argued in |5^, that the mixing IR divergences are cancelled in gauge invariant quantities. It 



is shown in |2^, that the mixing divergences are also cancelled completely in the old-fashioned 



Hamiltonian theory for the coupling constant renormalization. In x"'"-ordered Hamiltonian the- 
ory, the mixing divergences cannot and should not be removed as indicated by the negativity 
of the second order correction of the wave function renormalization constant, that must result 
in a physical theory |5^]. 



Zhang and Harindranath p4[ performed the similar calculations of self-energy mass oper- 



ators in light-front QCD, using perturbation theory frame. They have tested different types 
of regulators to regularize the transverse UV divergences; the regulator for IR divergences 
was chosen in accordance with the prescription for boundary terms, arising in the light-front 
Hamiltonian (see chapter |^). In the scheme discussed above the transverse regulator arise au- 
tomatically at the stage (2), where the unitary transformation, dictated by flow equations, is 
performed. This is important, since operator classiflcation and renormalization itself are driven 
by the scaling in transverse direction. The problem, that still remains, is the violation of rota- 
tional and gauge symmetries by this regulator, in other words, by the flow equations. We leave 
this problem for future investigation. 

Further the detailed calculations of electron and photon self energies and physical masses, 
using the program outlined above, are presented. 

7.2 Flow equations for renormalization issues 

As was discussed above and in chapters 3 and 4, the commutator [r^'^^^ Hee-^ contributes to the 
self-energy term, giving rise to the renormalization of fermion and photon masses to the second 
order. The flow equation for the electron (photon) light-cone energy p~ is 

J, ^ ['/ ; -'^ee^J ^ self energy i V ' •"/ 

where the matrix element is calculated between the single electron (photon) states < p', s'|...|p, s >. 
We drop the flnite part and deflne bp^ = p^{lx)— < \Ho\ >■ Integration over the flnite range 
gives 

X - A - /■'' / r (1) H- 1 ^ Ai' (^Sa(p)-(5Sa(p)) 

^P\ -^Pa= < H ,Hee-y\ >self energy dl = , (7.9) 

JIa P 

that deflnes the cutoff dependent self energy 6Tjx{p)- The mass correction and wave function 



renormalization constant are given correspondingly, cf. [^ 



fil 



d5p 

1 + 



dp 



(7.10) 



The on-mass-shell condition is defined through the mass m in the free Hamiltonian Hq. 

We show further, that to the second order O(e^) the electron and photon masses and 
corresponding wave function renormalization constants in the renormalized Hamiltonian vary 
in accordance with the result of 1-loop renormalization group equations. This can serve as 
evidence for the equivalence of the fiow equations and Wilson's renormalization. Therefore 
we have rewritten the mass correction 6m\ through the self energy term, arising in 1-loop 
calculations of ordinary perturbative theory. The negative overall sign stems from our definition 
of the fiow parameter, namely for A/ > we are lowering the cutoff dl = —-^dX. 

We start with the bare cutoff mass m\ = rn^ + 5Mp^ , where 5M)^ is the second order mass 



counterterm. According to eq. ( [7.9| ),eq. ( |7.10| ) the electron (photon) mass runs 



ml = ml- [(5S,(m2) - 5T.^{m^)] (7.11) 

defining, due to renormalizability, the counterterm 5M)^ = 5m\ = —STi^^m^) and the depen- 
dence of the renormalized mass on the cutoff A 



m 



l = m^ + 5ml = m^-5Ex. (7.12) 



We calculate explicitly the self-energy term. The electron energy correction contains several 
terms 

6p~, =< p', s'\H - Ho\p, s >= (t^ ^Pxn) ■ S^'\P - P)Sss' . (7.13) 

The first term is induced by the fiow equation in single electron sector, namely comes from the 
commutator [ri{l),Hee'y] 



^PlX = - < [V^^\ Hee^] >self energy dl' = ^-^ ; (7.14) 

Jlx P 

it reads, cf. eq. (|C.10| ) in Appendix 0, 



2 / d'^k^dk+e{k+)^ 



xr'(p - k,p, -k)r(p,p- k, k) ^.-\ ^.y X (-fl) ■ (7.15) 

This term explicitly depends on the cutoff A (/ = 1/A^) through the similarity function, that 
plays the role of a regulator in the loop integration 

i?A = &,. = exp|-2(^)'|. (7.16) 

Eq. ( [7.15|) corresponds to the first diagram in fig. (9). 

Two instantaneous diagrams, the second and third in fig. (9), contribute the cutoff indepen- 
dent (constant) terms. They arise from normal-ordered instantaneous interactions in the single 
electron sector and can be written as 

Spnx = Spnil = 0) = c„ < do-' > Vr\l = 0) (7.17) 

fi9 



where n = 2, 3 corresponds to the second and third diagrams in fig. (^, c„ is the symme- 
try factor, < 00^ > stands for the boson (n = 2) and fermion (n = 3) contractions (i.e. 
< (ika^ >= 6{k^)/k^ and < 6p6j" >= 6{p~^)), and l^"**(/ = 0) arises from normal-ordering of 
if 6677 for n = 2 and of -ffeeee for n = 3 (eqs. ( p.29| ) and (|3.3CI|) ). 

These two diagrams 6pn{l = 0) define together with the first one 6pi{l = 0) the initial 
condition for the total energy correction, eq. ( [7.13 ). 

Since the diagrams n = 2, 3 come from the normal-ordering canonical Hamiltonian at / = 0, 
they must accompany the first diagram for any fiow parameter I. In what follows we use for 
the instantaneous terms the same regulator R, eq. ( [T.KJI ) 



SP2X 



(Pk^dk+ e{k+) a'a' 
2(27r)3 k+ |j^ 



Q I (X K (IK /^/,_L\-L 

^(2^'^') 2 



k+] 
1 



x(-i?) 



[p+ - k+f [p+ + k+Y 



:-R) 



(7.18) 



We define the set of coordinates 

k+ 
p+ 



X 



k 



(7.19) 



where p = (p^,/?"*") is the external electron momentum. Then the electron self energy diagrams, 
fig. (9), eq. ( |C.14| ) in Appendix y, contribute 
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(7.20) 



and 
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(7.21) 



for details we refer to Appendix y. Note, that the transformation in the integrals over x is 
performed before the regulator is taken into account [^]. (In the second integral the electron 
momentum is replaced by the gluon one due to momentum conservation). The brackets '[ ]' 
denote the principle value prescription, defined further in eq. (|7.27| ). 

The loop integral over k eqs. ( [7.20| ) and ( |7.21| ) contains two types of divergences: UV in 
the transversal coordinate k-^ and IR in the longitudinal component k~^. The physical value 
of mass must be IR-finite. We show, that the three relevant diagrams together in fact give an 
IR-finite value for the renormalized mass; this enables to determine counterterms independent 



F,?, 



of longitudinal momentum. In the wave function renormalization constant, however, the IR- 
singularity is still present. 
Define 

^1 = ^, (7.22) 

where P = {P^, -P"*") is the positronium momentum, p the electron momentum. The transversal 
UV divergence is regularized through the unitary transformation done, i.e. by the regulator R, 
eq. (^) 




^A = exp - -^ ^ eiX% - |Ap,fc|) , (7.23) 



where the cutoff is rescaled and defined in units of the positronium momentum P"*", namely 
A -^ v^A^/P+, and Ap^k = p^ — k^ — {p — k)^ = Ap^k/p^- The rude approximation for the 
exponential through a 6'-function changes the numerical coefficient within a few percent; nev- 
ertheless it is useful to estimate the integrals in eqs. (|7.20| ) and (|7.21|) in this way analytically. 
From eq. ( |7.23| ) we have for the sum of intermediate (electron and photon) state momenta (the 
external electron is on-mass- shell p^ = m?) 



^±2 ^±2 _|_ ^2 



[x] [1 — x] 

giving for the regulator 



< \^5i + m^ (7.24) 



R\ — (^{f^Xmax ^ ) ^{^XmaxJ 

41. = xil-x)X%-x'm' (7.25) 

and ^(fi^Amax) leads to the additional condition for the longitudinal momentum 

U j^ X 2i •'^max 

1 + m^/[ydi) 

implying that the singularity of the photon longitudinal momentum for x — >■ 1 is regularized 
by the function Rx- This is the case due to the nonzero fermion mass present in eq. ( [7.24| ) for 



the intermediate state with (1 — x) longitudinal momentum. The IR-singularity when a; — >■ 
is still present; it is treated by the principle value prescription P^ 



^ ^ ' ^ + 7^^-^ ) ' (7-27) 



[k+] 2\k+ + ieP+ k+-ieP+ 

where e = 0+, and P~^ is the longitudinal part of the positronium momentum (used here as 
typical momentum in the problem). This defines the bracket '[ ]' in eqs. ( |7.20|) and ( [7.21| ) 



[x]-2[x + ^f^^x-^fJ ' ^^''^^ 
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Making use of both regularizations for transversal and longitudinal components, we have for 
the first diagram, eq. (|7.20|) , 
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Note, that the third term has the mixing UV and IR divergences. Combining the three relevant 
diagrams, fig. (9), and integrating with the common regulator, one obtains for the electron 
mass correction 



5m\ 
5m\ 



p^{5pi + 5p2 + (5p3)|p2=^2 = 
"2 r /\^5, + m^' 
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The mass correction is IR-finite (that gives rise to IR-finite counterterms) and contains only a 
logarithmic UV- divergence. Namely, when \5i ^ A ^ m 



Smi 



3e2 2, A' 
-m In — - . 



87r2 



m^ 



(7.31) 



It is remarkable that we reproduce with the cutoff condition of eq. ( [7.24[ ) the standard result 
of covariant perturbative theory calculations including its global factor 3/8. As was mentioned 
above, the difference in sign, as compared with the 1-loop renormalization group result, comes 
from scaling down from high to low energies in the method of flow equations. 

The similar regularization for the intermediate state momenta in the self-energy integrals. 



called 'global cutoff scheme', was introduced by W. M. Zhang and A. Harindranath [^]. In our 
approach the UV-regularization, that defines the concrete form of the regulator R, arises nat- 
urally from the method of flow equations, namely from the unitary transformation performed, 
where the generator of the transformation is chosen as the commutator 77 = [Hd,Hr]. Note 
also, that the regulator R, eq. ( |7.23| ), in general is independent of the electron momentum p"*" 
(rescaled cutoff A^i — >• A), and therefore is boost invariant. 

For the wave function renormalization constant, eq. ( fT.lOl) , one has 



d6p 



dp 



■^TT^ Jo 



dti^. 



2i-2 + a; x{l - x)2m^ 



kI + f{x) (4 + f{x)y 



x(-i?) 



(7.32) 



that together with the regulator i?, eq. ( [7.23| ), results 
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As XSi — i> A ^ m the function F tends to a constant 
F\ -C- 2 + "' 

-r |A>m — O — — - + — . 

Therefore, by dropping the finite part, we obtain 
-2 '' A2 /3 . 1\ . 1 
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(7.35) 



where we have rescaled 4 — > e. The electron wave function renormahzation constant contains 
logarithmic UV and IR divergences mixed, together with pure logarithmic IR divergences. We 
mention, that the value of Z2 is not sensitive to the form of regulator applied; the same result 
for Z2 was obtained with another choice of regulator pi|. 

We proceed with renormahzation to the second order in the photon sector. The diagrams 
that contribute to the photon self energy are shown in fig. (10). The commutator [r]^^\ H^ 
corresponding to the first diagram, gives rise to (eq. ( p.22| ) in Appendix y) 
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where momenta are given in fig. (10), and the regulator is 



R\ 



fq±,\ = exp <^ -2 
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(7.37) 



In full analogy with the electron self energy this also defines the regulator for the second diagram 
with the instantaneous interaction, see fig. (10), 
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k= {{1- x)g+, (1 - x)q-^ + K-^) 
{q — k) = (xg"*", xq^ — k"*") , 



(7.39) 



where q = {q^,q^) is the external photon momentum. Then two diagrams contribute (for 
details see Appendix y, eq. ( |C26D ): 
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Note, that the transformation in the second integral is done before the regularization (by 
regulator the R) is performed [^ . 

Making use of the same approximation for the regulator as in the electron sector, we obtain 
for the sum of intermediate (two electron) state momenta 
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where the photon is put on mass-shell g^ = and the rescaled cutoff A — > a/2A^/P+ has been 
used. This condition means for the transversal integration 
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and for the longitudinal integration 
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where the approximate value is used when m <^ X. This shows that the condition of eq. ( 7.41| ) 
for two electrons with masses m removes the light-front infrared singularities from x —>■ and 
X — > 1. Thus, both UV and IR divergences are regularized by the regulator R, eq. ( [7.42| ). 
The mass correction arising from the first diagram, eq. ( [7.40| ), is 
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Combining together both diagrams with the same regulator, eq. (|7.40D, we obtain 
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where r is defined in eq. ( |7.43 ). The result shows that the mass correction involves the quadratic 
and logarithmic UV divergences, i.e. as \52 -^ K^ m, 
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The wave function renormalization constant is defined through 
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that, with the regulator R, eq. ( [7. 421 ), results 
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The photon wave function renormahzation constant contains only logarithmic UV divergence, 
indeed as X62 — > A ^ m 

e^ 2 A^ 
" - ' In^) (7.49) 



StT^ 3 TTl^ 

and is free of IR divergences (as is expected from the form of the regulator R, eq. ( [7.41| )) 



7.3 Mass renormalization 

Following light-cone rules the perturbative energy correction of the electron with momentum 
p, coming from the emission and absorption of a photon with momentum k, is 



^^^^ " J 2(271)3 k^^^P^ ~ ^^'>9p-k,p,x^\{P - k,p, -k)gp,p^k,xT\{p,p- k, k) 



dPk^dk+e{k+) 
2(27r) 

where 5'ee7-coupling constant restricts the energy of the photon. Making use of the explicit form 
for the coupling, one has 

2 f d^k^dk+ e{k+' 



xrl(p - k,p, -k) r\(p,p - k, t)— -pJ-^— -— X (R) , 

where R = f^j^y^ plays the role of regulator. This expression coincide up to the overall sign with 
the energy correction obtained in the previous section from the flow equations method. 

Two instantaneous diagrams, arising from the normal-ordering Hamiltonian, must be added 
to the first term with the same regulator R. Then the full perturbative energy correction 

5pl = Spix + SP2X + Sp3\ is 

5px = -hi (7.52) 

where (5p^ is defined in eq. ( |7.13| ). This means for the perturbative mass correction 



and the self-energy term 5Sa is given in eq. ( |7.3CI| ) 



We combine the mass operator, renormalized to the second order, eq. (|7.12|) , and the per- 
turbation theory correction, eq. ( [7.53| ), to obtain the total physical mass to the order O(e^) 



ml = m{ + 6m' = {m' + 6Ex) - 6^x = m' + O(e^) . (7.54) 

This means, that to the second order O(e^) the physical electron mass is, up to a finite part, 
equal to the electron mass, that stands in the free (canonical) Hamiltonian Hq. 

Along the same line one can proceed for the photon mass. One finds in the second order 
in coupling, that the photon mass counterterm, obtained from the flow equations eq. ( [7.44]) - 



eq. ( 7.46| ), is equal up to the overall sign to the perturbative photon mass correction. This 



means, that in the order O(e^) the physical photon mass is equal to the photon mass term in 
the free Hamiltonian, i.e. 

ml, = + Oie') (7.55) 



here the photon mass in the canonical QED Hamiltonian is equal to zero to preserve gauge 
invariance. 

At the end we note, that the similarity function fp^pfX, restricting the electron-photon vertex, 
plays the role of UV (and partially IR) regulator in the self energy integrals. This means, that 
the regularization prescription of divergent integrals follows from the method of flow equations 
itself. Moreover, the energy correction (i.e. mass correction and wave function renormalization 
constant), obtained from the flow equations, coincide up to the overall sign with the 1-loop 
renormalization group result. This is the remarkable result, indicating to the equivalence of 
flow equations and Wilson's renormalization. 
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Chapter 8 

Conclusions and outlook 



In this work we applied the method of flow equations to QED on the light-front. We have 
outlined a strategy to derive an effective low-energy Hamiltonian on the light-front by means 
of flow equations, considered in the light-front dynamics. Application of the flow equations 
with the condition, that particle number conserving terms are considered diagonal and those 
changing the particle number off-diagonal led as in other cases to a useful effective Hamiltonian. 

The main advantage of this procedure as compared with the similarity renormalization 
of Glazek and Wilson ^ is, that finally states of different particle number are completely 
decoupled, since the particle number violating contributions are eliminated down to A = 0. Thus 
one is able to truncate the Fock space and the positronium problem reduces to a two particle 
problem, which was analyzed analytically (since in leading order one obtains the nonrelativistic 
Coulomb problem) in the chapter 4, and numerically in the chapter 5 for positronium bound 
states. 

The effective Hamiltonian, obtained by the similarity transformation, is band-diagonal in 
the energy space. The width of the band A introduces the artificial parameter in the procedure, 
which is defined from the physical reasoning (A is low enough to neglect the contribution of high 
Fock states, but is restricted from below to stay in perturbation theory region). Flow equations 
as used here with the particle number conserving part of Hamiltonian to be diagonal, have no 
additional parameter and converge well as A ^ to the effective Hamiltonian, which is 
block-diagonal in particle number and are used therefore directly for the numerical calculations 
of the spectrum. 

The procedure of elimination of nondiagonal blocks, that change the number of quasiparti- 
cles, is performed not just in one step as in the method of Tamm-Dancoff truncation but rather 
continuously for the states with different energies in sequence. This is the main advantage of 
the proposed method as compared with Tamm-Dancoff truncation, the possibility to perform 
simultaneously the ultraviolet renormalization of the initial Hamiltonian. In general, in the 
definite order of perturbation theory all counterterms, associated with canonical operators of 
the theory and also with possible new operators induced by unitary transformation, can be ob- 
tained in the procedure |^3[. Since different sectors of the effective Hamiltonian are decoupled, 
one does not encounter the usual difficulties of Tamm-Dancoff truncation and the methods 
related to it. Namely, the counterterms to be introduced are 'sector-' and 'state-' independent 

If one goes beyond the tree approximation then one obtains terms with ultraviolet di- 
vergences which have to be renormalized. In the second order in coupling the electron and 
photon divergent mass corrections are generated by flow equations. There are the same type 
of UV-divergences as obtained in the covariant perturbation theory, the IR-divergences in the 
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longitudinal direction and mixed UV- and IR-divergences in the mass terms. The two lat- 
ter types of divergences are specific for the light-front QED calculations. They signal, that the 
boost invariance is broken, since an explicit dependence on the longitudinal IR-cutoff occur. To 
preserve boost invariance in the renormalized Hamiltonian we take into account the diagrams 
arising from the normal ordering of instantaneous interactions. Using then fiow equations and 
coupling coherence we obtain the counterterms for electron and photon masses, that are free 
from IR-divergences. 

Simultaneously also terms describing interactions between more than two particles are gen- 
erated. In this approach we were not faced with infrared problems, except for the longitudinal 
IR-divergences, which arise due to the light-front gauge formulation and must be treated prop- 
erly. These divergences, arising on the level of effective QED Hamiltonian, show that some 
symmetries of the initial theory, that are not manifest on the light-front, are broken. The first 
problem in the light-front field theory is that whenever the generator of a symmetry is dynam- 
ical (contains interactions) it is practically impossible to monitor and maintain that symmetry 
at each step of a calculations- unless of course one can solve the theory exactly. This concerns 
parity and rotational invariance, that are not manifest on the light-front. The second problem 
is that whenever we use the Hamiltonian technique, the regularization by introducing bare cut- 
offs breaks the gauge invariance (and also Lorenz covariance) , and forces the bare Hamiltonian 
to contain a larger than normal suite of counterterms to enable a finite limit as the cutoffs are 
removed. One way is, that the counterterms are then adjusted to reproduce physical observ- 
ables and to restore the symmetries broken by the cutoffs. The other way is to find the gauge 
invariant procedure for regularization of Hamiltonians. The attempt in the latter direction was 
made by Brodsky, Hiller, McCartor |^. Solving the Yukawa theory they have tried to preserve 
more symmetries by using Pauli-Villars procedure for regularization. 

In our calculations, the effective electron-positron interaction, obtained in the second order 
in coupling, contains the IR-divergent term, describing instant emission and absorption of 'lon- 
gitudinal' photon. The physical reason for its appearance is the violation of Lorenz and gauge 
symmetries by the derivation of the effective, renormalized Hamiltonian. We use the symmetric 
cutoff condition for IR-divergences in the effective interaction to cancel their contribution to 
the mass of positronium. We find, that in this case the rotational symmetry is restored on the 
level of positronium mass spectrum. 

In order to solve the fiow equations analytically we were forced to apply in this work the 
perturbation theory expansion. One is able to improve this approach systematically by going 
to higher orders in the coupling. It is a remarkable result, that the effective electron-positron 
Hamiltonian, obtained in the second order in coupling, gives the correct Bohr spectrum and 
hyperfine splitting for positronium. 

We consider fiow equations as a method which can also be used beyond perturbation theory 
in a self-consistent way. Examples in solid-state physics are the fiow of the tunneling-frequency 
in the spin-boson model [^] and of the phonon energies in the electron-phonon coupling |T0| ]. 
Due to the fiow the couplings decay even at resonance. 
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Appendix A 

Calculation of the commutator 

[t]^ \l)^Heej] in the electron-positron 
sector 



Here we calculate the commutator [q^^'^l), Hee-y] in the electron-positron sector. The leading 
order generator i]^^^ is: 

^^^HO = E J d^Pid^P3d^q{v*p,p,i^yxag + np,p,{l)e'^atg) (b+^bp^ + b+Jtp^ + d-pAi + d-vJ-p,) 

xxirKpi.P3,-g)xsi'^9,-(pi-P3) > (A-1) 



i,S3 

where 



J- ^9piP'j, 



VpiPsK^) ~ ^PlP3 ' 9piP3 ~ \ ' rll '' v^-^J 

^piP3 = Pi ~ Ps ~ (Pi ^ P^)i and the electron-photon coupling 

Hee^ = J2 J d^P2d^Pid^qig*p^p^{l)e{,ag> + gp2pSyx*dtg>) (6+ 6p2 + b^Jtp^ + d-pj>p2 + d^pjtpj 
AS2S4 

^xtJliP2,P4,-q')Xs2Sq',~{P2-P4) > (A.3) 



where 

q+ pj pj 

and the tilde- fields are defined in eq. (|3.31|) . The mass term nip^p^il) starts to depend on / in the 
second order in coupling constant. We omit the dependence oinip^p^ and hence the dependence 
of vertex matrix element T^{pi,P2iq) on / in calculation of commutator [ri^^\l),Hi,e^] to the 
leading (second) order. Further we use the identities for the polarisation vectors and spinors 

^£r4 = ^^^ xtxs' = 5ss'. (A.5) 

A 



Using the commutation relations, eq. ( p.23|) , and identities eq. ( |A.5| ) we have 



h'^Hl) H ] = -i-v q* (^^Pt-Pt) . 0{pt-pt) \ .j^Q. 

[I I yi),i2ee-y\ „ '/piP3yp2P4 ^+ ^+ ^'ipiP3yp2P4, + + y-^-^) 

-^ \ Pi ~ Ps P3 ~ Pi / 

^. [^ Up,^U,_p^U,^p^Up^ Up^U_p_^U^p^Up2 T Up,^U,_p^Li^p^Up^ -r Up^U,_p^U,^p.j^Up^) . 

X(xir;(pi,P3,Pl -p3)Xsi) (X^4r;(p2,P4,P2 -P4)XS2) (5pi+p2,p3+p4 , 

7r>. 



where the first two terms of the field operators contribute to the exchange channel, and the 
next two to the annihilation channel. We take into account both s- and t-channel terms to 
calculate the bound states. The : : stand for the normal ordering of the fermion operators and 
(|) is the symmetry factor. The sum over helicities Si and the 3-dimensional integration over 
momenta pi, i = 1, ..4, according to eq. ( |3.32| ) is implied. We rewrite for both channels 

_p+) \'ipi,P'j.y-pi-p2 '/-p4,-p2ypi'P3j 



hH, 



+ 



29r 



/pi .p.^y- 



-(pt^Ps) ^''^P'l,-P2ypi,p3 'lpi,P3^-P4,-P2 

S h+ rl+ rl - h 

-'Pl +P2 ,P3+P4 "p3 S3 "'P4S4 "'P2S2 "Pl Si 



ee'y] 



(A.7) 



-M. 



(an)^l^^ J eip++p+) 
2ij 



(2) 1 (P++P+) 



in. 



'Pl,-P2fi'-p4,p3 



S{-{pt+Pt)) . 

-(Pt+Pt) 



~l /„+^„+\ V/-P4,P3ypi,-p2 



n-pi,p-i9pi-p2) 

~ Vpi-p29-p4,P3y 



where 



M 



M 



(ex) 
2ij 

(an) 
2ij 



vA*JA h+ rl+ rl - h 

/\U Up-^_(_p2,p3+p4 *^p3S3"'p4S4"'P2S2*^piSl 



(xir'(Pl,P3,Pl -P3)Xsi) (xJ,r-'(-P4, -P2, -(Pl -P3))Xs-4) 



(xir'(-p4,P3,-(pi +P2))x54) (x^2r-'(pi,-p2,Pi + P2)xsi; 



(A.8) 



The first term in the exchange channel with pi > p^ corresponds to the light-front time 
ordering xf < xf with the intermediate state Pf7 = p^ + {pi — Ps)^ + P2 ^ the second term 
Pi < Pt ^^^ ^1 > ^3^ has the intermediate state P/T = Pi — {pi — Ps)' + Pa- Both terms can 
be viewed as the retarded photon exchange. The same does hold for the annihilation channel. 

Consider only real couplings and take into account the symmetry 



V- 



P4.-P2 



-Vp 



'4,P2 1 



9- 



P4.-P2 



9p. 



'4,P2 



(A.9) 



Then ^353,^434! [r]*^^^ -^^667] \piSi,P2S2) , the matrix element of the commutator between the free 
states of positronium in the exchange and annihilation channel, reads 



<[ri^'\H, 



667 J 



> = 



^'^2m („+_„+) \Vpi,ps9p4,P2 ' Vp4,p29pi,: 



■P3) 



'2m 



(Pt+Pt) 



29p 



P2HP4-P3 



+ Vp4,-P39pi,-p2J 



(A.IO) 



where the conservation of '+' and ' ±' components of the total momentum is implied, i.e. 
Pi ~^pt = pt ~^Pt Siud pI+P2 = pi+pi- We rewrite this expression through the corresponding 
/-functions 



Vpi,P39p4,P2 ' Vp4,P2 9pi 



P3 



'^/pi.psvJ 



_^pi,P3 



Vpi~P2 

with A 



9p4-P3 "T Vp4-P39pi,-P2 
P2) 



dl 
"/pi,-p2l'j 



/P4,P2 VV ~r 



"/p4,P2l'J 



A 



P4,P2 



dl 



/pi.Pslv 



(A.ll^ 



A. 



dl 



/P4,-P3lv + 



"/p4,-P3lv 



A. 



dl 



Jpi,-P2\'') 



P1-P2 ""- ^P4,-P3 

This form in terms of the /-function is universal for all 



Pl,P2 = Pi -P2 - (Pl 

unitary transformations. 

We calculate the matrix elements M2ii, eq. (188), for both channels. Here we follow the 
notations introduced in p9 |. 



7."^ 



We make use of the identities 
with s = —s and xtXs' = ^ss'] also of 



xt^'Xs 



-V2. 



se. 



xt^^'x-s 



-V2. 



se: 



(A.12) 
(A.13) 



with e* = —Es and 5*5*, = —6g->. 

We use the standard hght-front frame, fig. (3), 

pi = {xP+, xP^ + /t^) , p2 = ((1 - x)P+, (1 - x)P^ - /t^) , 

P3 = {X'P+, X'P^ + k'^) , P4 = ((1 - X')P+, (1 - X')P^ - k'^) , 

where P = {P^,P^) is the positronium momentum. 

Then, to calculate the matrix element M2ii in the exchange channel, we find 



(A.14) 



P^[xtJ'iPuP3,Pl -P3)Xsi] 



Xs-i 



^ (k,± — K,', y Cr ■ k'i ,• ,• O" ■ /t_L . X — x' 

2^ ^ —a' + a' + im a 



{x — x') x' 



X 



XX' 



T^6s,s, + im —{-V2)siei6s, 

XX 



S3 5 



Xsi 

(A.15) 



and 



P'^[xt2^'{-P4, -P2, -{P4 - P2))X. 



Xs2 



S4j 



2 ■ h a + a 






im- 



X — x' \1 — X 1 — x' ) (1 — x)(l — x') 



-a 






T16.2S. + ^^ n _ ,.n _ ,,. i-^^>^^\A 



Xsi 

(A.16) 



where we have introduced 



Tl 



T' 



[K±-K\y , kUs2) , k'\(s2) 



X — x' 



:i -x) (1- x' 



(A.17) 



(K±-K\y kUsi) k'\(si] 



X — x' 



X 



x' 



and 



tx, \ \S ] — Kj I ~r '2' 5 S'ij t\j I 



(A.18) 



Finaly we result 



+2 i\Aex) 



P+^M. 



2ii 



{"siS30s2S4-'l ■-'2 dsiS2^SiS3ds2S42m, 

xx'{l — X){1 — x' 



(A.19) 



+im^/2{x' — X 
Whereas in the annihilation channel we calculate 

P^[xt'^'{-P4,P3,-{Pl+P2))Xs4] = Xt 



e _ r ^Lt^ c^ -I- a a — T-"- r-^ I 



-^a' + a'- 4 + ^m- 

X 



1 — x' x'{l — x') 



a 



■'S3S4 



XSA 

(A.20) 
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and 



X. 



S2 



a — a tm^. ra 

I — X X 

1 



TlSs-,s2 - im 



x{l — x) 



x(l — x) 
-V2)siel^6s,s2 , 



Xsi 



where we have introduced 



7^3 


= 


x' 


l-x' 


T' 




/.Ksi) 


K\{s^) 



1 — X 



X 



We finally have 



+2 i\Aan) 



P^^M. 



2m 



XX [I — x)[\ — x' 



+imy2 



Sl 



S3 



O.3S4OS1S, .^_ ,^3 ^si ''-3S4^sis.^,.^_^,.J4 ^S4 



(A.21: 



(A.22) 



(A.23) 
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Appendix B 

Matrix elements of the effective 
interaction. Exchange channel 



In this Appendix we follow the scheme of the work PO] to calculate the matrix elements of 



the effective interaction in the exchange channel. Here, we list the general, angle-dependent 
matrix elements defining the effective interaction in the exchange channel eq. ( |6.55D (part 
I), and the corresponding matrix elements of the effective interaction for arbitrary J^, after 
integrating out the angles eq. ( |6.6|) (part II). The whole is given for the similarity function 
/a(A) = ux{A) = 6{X^ — |A|) eq. ( p.34| ) with the sharp cutoff. The effective interaction 
generated by the unitary transformation in the exchange channel reads eq. (|6.58 ),eq. ( |6.61 ) 



reff _ _^2 . ^nu ^ ^. / t7(ai - ttg) , e{a2 - ttl) 



VIU = -e^ < 7^^ >9,A ' , ^' + 



Ai A 



2 



=2 ^ o/Mo/*^ 



1 f {ai - a2)9{ai - a2) (02-01)^(02-0-1' 



-e < ri > V,V.^2 y ^^ 



Ai A 



2 



-e^<rr>v.^i<..-.i(^^ + ^^) (B.i; 



where fig. (4) 



< 7^^ > U, = ^^P^^)^Ap^X) v{P2X)^MP2^.^p^2 (b 2) 

yPi yp'i' yp'2' 

q = p'l — Pi is the momentum transfer. One has in eq. 

Ai = oi — 2k±kj_ cos{{p — (p ) 

A2 = 02 — 2k±kj_ cos{{p — (p ) 

A = o — 2k±kj_ cos((y9 — (f ) 

k± = k±{cosip,smip) (B.3) 
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1- 




~ ^ 1 '2 , 2 ( 


X — X 

x)(l 


-X') 










kl 


+ k!l + (x 


-x') 
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— X 
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-xV 
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'"V 


(x-x')2 
-x)(l- 


X') 


kl 


, /n (X 

+ kl + ^ 


-x') 
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(a:1( 
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4' 


r '2 




1 

-x' 


iO 




1-x 
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{x - x'f 1 
2 [ 


1 
xx' 


(1- 


1 
x)(l 














-(ai + os) 





















(B.4) 

The energy denominator (A and a corresponding) in the case of perturbative theory is given 
for completeness. 

It is useful to display the matrix elements of the effective interaction in the form of ta- 
bles. The matrix elements depend on the one hand on the momenta of the electron and 
positron, respectively, and on the other hand on their helicities before and after the inter- 
action. The dependence on the helicities occur during the calculation of these functions 
E{x, k±; \i, \2\x' , kl; \[, \2) in part I and G(x, /c_l; Ai, A2|x', A;^; A'^, A2) in part II as different 
Kronecker deltas jESl. These functions are displayed in the form of helicity tables. We use the 



following notation for the elements of the tables 

Fi(l,2) ^ Ei{x,kr,x',kl); Gi{x,kr,x',kl) (B.5) 

Also we have used in both cases for the permutation of particle and anti-particle 

F;{x, k±; x', k'^) = Fs{l - x, -fc^; 1 - x', -kl) (B.6) 

one has the corresponding for the elements of arbitrary J^; in the case when the function addi- 
tionally depends on the component of the total angular momentum Jz = n we have introduced 

F,(n) = F,{-n) (B.7) 

B.l The general helicity table. 

To calculate the matrix elements of the effective interaction in the exchange channel we use 
the matrix elements of the Dirac spinors listed in Table B.l [|^. Also the following holds 
vx'iph^'vxiq) = ux{qh"'ux'ip). 



We introduce for the matrix elements entering in the effective interaction eq. ( B.l 



2E^^\x,kr,Xi,W,kl;X[,Xl) = <j^'Y>9,.u = 

= -< 7+7- > +2 < 7^7+ > - < 7? > - < 72 > 
2E^'\x,kr,X,,X2\x',kl;X[,X',) = < YY > V.V.^^ =< 1^1^ > ^ (B.8) 



where 



^ ^^ u ^_ u{x, kr, Ai) Y u{x', kl; X[) v{l - x' , -kl; A^) Y ^(1 - x, -k±; A2) .^ ^. 

Jxx'{l — x)(l — x') 
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M 


^uik',X')Muik,X) 


7+ 


24 


7" 


2 


"(m^ + A;iA;le+*^(^-^')) ^ - mA (A;le+*^'^' - k^e+'^^) 6^y 




A;+A;'+ 


7i 




7i 


( h' , k , \ /I 1 \ 



Table B.l: Matrix elements of the Dirac spinors. 



final : initial 


(A'i,A'2)=TT 


(A'i,A'2)=Ti 


(A'i,A'2)=iT 


{K,K)=ii 


(Ai,A2)=TT 


i^i(l,2) 


E*,{1,2) 


^^3(1, 2) 





(Ai,A2)=Ti 


E*{2,1) 


i^2(l,2) 


^4(1,2) 


-^3(2,1) 


(Ai,A2)=iT 


^3(2,1) 


^4(1,2) 


^^2(1, 2) 


-^3* (2,1) 


(Ai,A2)=ii 





-^3(1, 2) 


-E*{1,2) 


i?i(l,2) 



Table B.2: General helicity table defining the effective interaction in the exchange channel. 

These functions are displayed in the Table B.2. 

The matrix elements E^ {1, 2) = E^ {x, k±;x',k'j_) (n = 1, 2) are the following 



e[^\x, k±; x', k'^) = m^ \ — ; + 



+ 



k\k ^ 



A{ip~ip ) 



xx' (1 — x)(l — x')y xx' {\ — x^{\ — x'^ 



E^^\x,kr,x',k'j 



m^ I — - + 



+ A;^ ^T^ + k7 ^ 



xx' (I — x)(l — x') J x{l — x) x'(l—x') 



+ kxk 



Jif-v ) 



-i{^-ip ) 



n(l) 



xx' (1 — a;)(l — x') 

m I , , i,J . 1 — x' 



E^'>{x,kr,x',k^) = [k^e''^ -k^- 

xx' \ 1 



lip 



X 



J-^A ['^ J '^_L 1 '^ •) '^ \_ 



-m 



{x — x') 



,/\2 



xx'(l — a;)(l — x') 



(B.IO) 
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and 



El (x, k±; x\ kj_) = E2 (x, k±_] x , kj 



(2), 



{x — x') 



,/^2 



-£^3 (x, k±; x\ kj_) = E[ (x, k±; x', kj_) = 



(2), 



(B.ll^ 



B.2 The helicity table for arbitrary J^. 



Following the description given in the main text eq. ( p. 61) we integrate out the angles in the 
effective interaction in the exchange channel eqs. ( |6.61|) and ([B.l| ). For the matrix elements of 
the effective interaction for an arbitrary Jz = n with n G Z G'(x, k±\ Ai, A2|x', fc^; A'^, A2) =< 
x.,k^]Jz.i\i.i\2\VlJ\x\k'^]J'^,\'^,\'2 > \exch in the exchange channel one obtains the helicity 
Table B.3. 



final : initial 


(a;,a'2)=tt 


(A'i,A'2)=n 


(A;,A'2)=iT 


(A;,A'2)=ii 


(Ai,A2)=TT 


Gi(l,2) 


G^(l,2) 


G3(l,2) 





(Ai,A2)=n 


G*(2,l) 


G2(l,2) 


^4(1, 2) 


-^3(2,1) 


(Ai,A2)=iT 


G3(2,l) 


^4(1, 2) 


^2(1, 2) 


-Gl{2,l) 


(Ai,A2)=ii 





-^3(1, 2) 


-G*(l,2) 


Gi(l,2) 



Table B.3: Helicity table of the effective interaction for J^ = ±ra, x > x'. 



Here, the functions Gj(l, 2) = Gi{x, k±; x', k'j_) are given 

'1 1 ~ 



Gi{x,k±;x',kj_) = m 



, /r2ta,a,(|l — nl 

XX' (l-x)(l-x')/ 



-L*^! T /I i\ K^l ^^2 T /I i\ 

+ xx'(l - x)(l - ^0'"'--^'"'^ + (^3^'"'--^"' - "") 
G2(x,fc,;x',fcl) = {m- (^ + (, _ ^)', _ ^.) ) + ^1^(1^ + ^^'^7(Y^)^^^^.^.(NI) 

yxx (1 — x)(l — X ) J 

(x — X Y 



G^{x, k\ \x' ,kt 



., ,VJ_, ^ , i 



m 

XX 



- I fc^/nta, 02(11 + n\) - ^±Y^^ — Inta^a2{\n\) 



'(l-x)(l-x') 
where we have introduced the functions 



G'4(x,/cx;x',/cl) = -"^' 3.^/n - a;Vi _ ^// ^^"i^^d^D 



(B.12) 



Inta^aiin) = 0{ai - a2)Inta^{n) + 61(02 - ai)Inta2{n) 



IntaXn) = -i-A{ai)) 



a , . , .._„^i ( B{ai) 



k\ k\ 



7P 



A{ai 



al - Aklk'l 



B{ai) = -{l-aiA{ai)) 



(B.13) 



and the functions Oj, i = 1, 2 are given in eq. (|B.4|) . 



The following integrals were used by the calculation of the matrix elements [|30 

cos{n{ip — ip')) 



— / d(f dip' - 
27r Jo Jo tti — 2k±kj_ cos{ip — Lp') 
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(B.14) 



The condition on the parameter space (x, kj_) due to the 6'-function, namely 6{ai — 02) (i.e. 
when ai > 02) reads 



{x - x') (x(l - x){k'l + m^) - x'{l - x'){kl + m^)) > 



Making use of the coordinate change eq. (|6.71 ) this is equivalent to 
fi cos 9 /i' cos 9' 



\fWTm^ J ^'2 + rn 



{jj, - jj,') <0 



(B.15) 
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Appendix C 

Fermion and photon self energy terms 



We calculate here the fermion and photon self energy terms, arising from the second order 
commutator [r]^^\ Heey]- 

I. We first derive the electron self energy terms. Making use of the expressions for the 
generator of the unitary transformation r]^^^ defined in eq. ( |4.18| ) and of Hee-y, eq. ( p.27| ), we 



obtain the following expression for the commutator in the electron self energy sector 



n\VpiP29p2Pl Vp2Pl9piP2) 



Pi -Pi 
where 



oiPt) ^J ''p o{pt)b;^b,,xt2Xs2 
P2 ~ Pi 



M2ijipi,P2) = ^\Pl,P2,Pl - P2)PiP2,Pl,P2 - Pi) (C.2) 

and the momentum integration over pi,p2 is implied; 1/2 stands as the symmetry factor. The 
matrix element of the commutator between the free fermion states is 



< Pi, Sl\['n^'^\ Hee-y]\pi, Si >selfenergy 

e{pt-pt ) 

Ip2 "■- ■-' ■—■--- - Pi -Pt 



- / iVpiP29p2Pi - Vp2Pi9prP2) ^iPt) — T 1- M2iiipi,P2) , (C.3) 

JV2 jOl 09 



where the integration /^ is defined in eq. (|3.32|) . We use the expression for the generator rj 
through the coupling, namely 

^ I "'9p2Pl _|_ ^9piP2 \ /p A\ 

VpiP29p2Pl ~ Vp2Pl9piP2 ~ 'T I 9pip2 'Tj I" 9p2Pl 'Tj I • l^-^J 

Change of the variables according to 
Pi = P 

P2 = Pk 

P1-P2 = k (C.5) 

brings the integral in eq. ( |U.3| ) to the standard form of loop integration 

f 6(k^) 

- J {7]p^p-k9p-k,p - 'np-k,p9p,p-k) d{p'^ - k'^) ,^ M2ii{p,p- k) . (C.6) 

81 



According to eq. (|7.9|), the integral /^^ of the commutator [ri^^\ Hge-y] defines the difference 
between the energies (or energy corrections) Spl^^ — Sp^j^. Making use of 
/•'a ^ 1 

/ "' [VpiP29p2Pl Vp2Pl9piP2) ~^ ~ 7 \^ \9pi,P2,A9p2,Pl,>' 9pi,P2,^9p2,PlA)\^' ' ) 

"''a Pi — P2 — [Pi — P2) 

we have the following explicit expression: 



^PiA - ^P 



lA 



2(27r) 



k+ 



-e(p+ - k- 



p — k — {p ~ k)' 



(C. 



xV{p - k,p, -k)T\p,p- k, k) 



exp 



^p,p—k 

X 



exp 



A ^ 2~ 

^^p,p— A; 

A 



where the solution for the ee7-coupling constant was used. Therefore the electron energy 
correction corresponding to the first diagram, fig. (8), is 

SPIX = e^/^^S^^%^-^^) (C.9) 



2(27r)3 k+ 
xT\p - k,p, -k)T\p,p- k, k) 



1 



-R) 



p — k ~ {p — k) 
where we have introduced the regulator R, defining the cutoff condition (see main text), 



i? = exp <^ -2 






(CIO) 



(note that Ap ^ = App_fc). To perform the integration over k = (A;+, A;-*-) explicitly, choose the 
parametrization 

k+ 



X 



pi 



k = {xp^ , Xp^ + K^) , (C.ll) 

where p = {p~^,p~'~) is the external electron momentum. Then the terms occuring in 6pi)^ are 
rewritten in the form 

P(p _ k,p, -k)T\p,p- k, k) = (^+).(Vx)^ ((4^^ - 4i + 2) 4 + 2mV) 
Ap,,„, = p--k~~{p~ ky = -,:^^{x{l - x)p' - 4 - xm') = ^ . (C.12) 

Therefore the integral for the electron energy correction corresponding to the first diagram of 
fig. (8) takes the form 



p~^6p 



lA 



^TT^ JO 



dx / dn^i 



2^2 



;^- | + l)4 + m^x 
(l-x)(Ki + /(x)) 



X 



-R) 



iir'^ Jo 



dx dn, 



X 



2 2 

p — m 



4 + /(x) \[x] 



— -2 + X 



2w? 



4 + /(a;) VN' [1-a;]/. 



(C.13) 



x(-i?) 



where 



X 



(C.14) 
was introduced (see main 



/(x) = xm^ — x(l — x)p^ . 

In the last integral the principal value prescription for A as 

text), to regularize the IR divergencies present in the longitudinal direction. 

We thus have derived the expression for the energy correction which has been used in the 
main text. 
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II. We repeat the same procedure for the photon self energy. The second order commutator 
[ri^^\ Hee-y] gives the following expression in the photon self energy sector 

{Pl - P2 



r)\VpiP2 9p2Pl Vp2Pl9piP2) ' 



oiptwi-pt) '!:i 'i: <«-.4*4 (C.15) 



at +\ot +\"\P2 Pl) + i j* 
-0[-Pi)d{P2) -JZT I+T «g «9^ a4 



TrM2ij{pi,P2)S, 



g-{pi-P2) 5 



(pt - Pl 

where M2ij(pi,P2) is defined in eq. ( |C2| ) and the trace acts in spin space; the integration over 
the momenta q, pi and p2 is implied. The matrix element between the free photon states reads 

< g, A I [r/^^^ , Hee^] \ q, A > self energy Sij (C. 16) 

= / iVpiP29p2Pi - VP2Pl9p^p2) di-ptWipt) TrM2ijipi,P2) (^g.-bi-pa) , 

q "'P1.P2 
that can be rewritten after the change of coordinates according to 

Pl = -k 
P2 = -{k-q) 
P2-P1 = q (C.17) 

in the following way 

-^ / {Vk,k~q9k-q,k - Vk~q,k9k,k~q) 0{k^)e{q'^ - k^) TrM2ij{k, k - q) , (C.18) 

q Jk 

where the symmetry 

V-Pl,-P2 ~ ~Vpi,P2 

9-pi,-p2 ^ 9pi,P2 [^.iJ) 

has been used. The integration of the commutator over / in the flow equation gives rise to 

(^,u - S,:^S'' ^ ^e^ I '^Hk-)e(,* - k^) ^__,_':\'^_,)_ (C,20) 



X 



Tr (r(fc, k-q, q)P{k - g, k, -q) 



exp { -2 ( ^)'| - exp |-2 (^'^^ ' 



This means for the photon energy correction 

5qi,5'' = ^e' f^^^^e{k^)0{q^-k^) (0.21) 



q+ J 2(27r 



X Tr {r{k, k - q, q)T^{k - q, k, -q)) ^_ _ ^_ _ ^^ _ ^^_ x (-i?) 
where the regulator R 



i? = exp|-2(^)'| (C.22) 

has been introduced. Define the new set of coordinates 

T = X 

q+ 

k = [[I - x)q+ , {1 - x)q^ + K^) 
q-k = {xq^ , xq^ - K^) , (C.23) 



where q = (g"*", g-*") is the photon momentum. Then the terms present in Sq-^^^ are 
r{k, k - g, q)r{k - g, k, -g) = ,',_^,. ((^^-^ + l) '^^^ + ^ 



(g+)2x{l-a;)2 

g+3;(l— a;) ' g^ 



A.-.. = g" - A:- - (g - A;)" = -^^l±^ + ^ = ^ • (C.24) 



The integral for the photon energy correction corresponding to the first diagram of fig. (9) takes 
the form 



+ x 


e2 


f\ 


1 oqix-- 


87r2, 
e2 






87r2, 


V 



_2 (2:^-2 + ^)4 + ^, 



dx / dKi^ ^^^^^x i-R) (C.25) 

dx I dK\ \ ^ (2x^ - 2a; + l) + "^ + ( -2 + — r ) Ix (-i?) 



k\ + f{x) V y [i_a;] \^ [x][l-a;] 

with 

f[x) =m^ - g\(l - x) , (C.26) 

and the principal value prescription, denoted by '[ ]', introduced to regularize the IR divergen- 
cies. 

This is the form of the photon correction used in the main text. 
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Figure 1: Flow equations perform the block-diagonalization of the bare Hamiltonian of the 
canonical theory Hb{A) into a Hamiltonian consisting of blocks with equal number of particles. 
For a finite value of A the matrix elements of the 'particle number changing' sectors are squeezed 
into an energy band with roughly \Ei — Ej\ < A (left hand side picture) and are eliminated 
completely as A ^ (right hand side picture). 

Table 1: The effective light front QED Hamiltonian matrix up to second order in e in 
the Fock space representation. The matrix elements of the 'diagonal' (Fock state conserving) 
sectors are unrestricted in the energy differences; the 'rest' (Fock state changing) sectors are 
squeezed roughly in an energy band of width A. Black dots correspond to zero matrix elements 
in order O(e^). Instantaneous and disconnected diagrams are not included. 

Figure 2: The matrix elements of the effective Hamitonian, obtained by the flow equations 
in the second order in coupling, with corresponding diagrams. The diagrams 2 — 5 belong to 
the 'diagonal' sector; the 1, 6, 7 correspond to the 'rest' sector (the 6, 7 diagrams are drawn 
schematically, namely the corresponding momentum change must be performed to get the real 
'rest' diagrams, depicted in Table 1.) The photon momenta are a;"'"-ordered, from left to right. 
The similarity function is chosen fp-p^^x = exp{—Ap^pJ\'^), where Ap^pj,x = I]pl —T,pJ (the index 
'i' denotes initial and '/' final states) and Ap-^p^x = p]^ —p2 —{pi—p2)~ , p~ = {p\+m\)/p'^ . 

Figure 3 : The effective electron-positron interaction in the exchange channel; the diagrams 
correspond to the generated and the instantaneous interactions. 

Figure 4: Another choice of the coordinates in the effective electron-positron interaction, 
that are used in the light-front integral equation. 

Figure 5: Positronium spectrum for —3 < Jz < 3, a = 0.3 and A = 1. The annihilation 
channel is not included. For an easier identification of the spin-parity multiplets, the corre- 
sponding non-relativistic notation Lf is inserted. Masses are given in units of the electron 
mass. 

Figure 6: Stability of positronium spectrum for J^ = 0, without annihilation interaction. 
Eigenvalues Mf for a = 0.3 and A = 1 are plotted versus A^, the number of Gaussian points. 
Masses are in units of the electron mass. 

Figure 7: Deviation of corresponding eigenvalues for J^ = and J^ = I {a = 0.3 and 
A = 1) with growing number of integration points A^. The graphs shown AM^ = Mf^{Jz = 
0) — M^{Jz = 1) for the states 1^5*1 (triagles), ^Pi (squares), and fPi (circles). 

Figure 8: Electron self energy term includes the diagram arising from the commutator term 
[r]^'^\ Hee'y] in the electron sector, and also two diagrams, the second and the third on the figure 
5, coming from the normal ordering of instantaneous interactions Hg^-yy and -ffgeee, resp. 

Figure 9: Photon self energy term includes the diagram coming from the commutator 
[r]^^\ Heej] in the photon sector, and also the diagram from the normal ordering of the instan- 
taneous interaction Hee-y-y- 
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